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ABSTRACT 



Context, Millimetric observations have measured high degrees of molecular deuteration in several species seen around low-mass 
protostars. The Herschel Space Telescope, launched in 2009, is now providing new measures of the deuterium fractionation of water, 
the main constituent of interstellar ices. 

Aims. We aim at theoretically studying the formation and the deuteration of water, which is believed to be formed on interstellar grain 
surfaces in molecular clouds. 

Methods. We used our gas-grain astrochemical model GRAINOBLE, which considers the multilayer formation of interstellar ices. 
We varied several input parameters to study their impact on water deuteration. We included the treatment of ortho- and para-states of 
key species, including IT, which affects the deuterium fractionation of all molecules. The model also includes relevant laboratory and 
theoretical works on the water formation and deuteration on grain surfaces. In particular, we computed the transmission probabilities 
of surface reactions using the Eckart model, and we considered ice photodissociation following molecular dynamics simulations. 
Results. The use of a multilayer approach allowed us to study the influence of various parameters on the abundance and the deuteration 
of water. Deuteration of water is found to be very sensitive to the ortho-to-para ratio of IL and to the total density, but it also depends 
on the gas/grain temperatures and the visual extinction of the cloud. Since the deuteration is very sensitive to the physical conditions, 
the comparison with sub-millimetric observation towards the low-mass protostar IRAS 16293 allows us to suggest that water ice is 
formed together with COt in molecular clouds with limited density, whilst formaldehyde and methanol are mainly formed in a later 
phase, where the condensation becomes denser and colder. 

Key words. Astrochemistry, ISM: abundances, ISM: clouds, ISM: molecules, Molecular processes, Stars: formation 



1. Introduction 

Understanding the formation of water is crucial, not only 
because of its primordial importance for life on Earth, but 
also because it is thought to be one of the most abun- 
dant oxyg en-bearing spec i es an d also one of the main gas 
coolants dCeccarelli et ail 1 19961: iKaufman & Neufeldl [l996; 
van Dishoeck et al. 2011). Interstellar water is believed to be 



formed mainly via three main mechanisms: 1) cold gas- 
phase chemistry, starting from the ionization of H2 by cos- 
mic rays, eventually leading to H30 + via ion-neu tral reactions 
that then recombine wi th electrons to form H2O dBate si [19861: 
Hollenbach et al. I I2009I) ; 2) on the surface of interstellar dust 
particles, via the hydrogenation of accreted atomic and molec- 
ular oxygen occurring at cold temperature s (Tielen s & Hagenl 
ll982tlCuppen &Herbstll2Tj07tlMivauchietal.ll2008l) : 3) warm 
gas chemistry, initiated by a few endoth ermic reactions involv- 
ing H 2 in warm gas (T > 250 K) dCeccarelli et all 1 19961: 
IKaufman & Neufeldfl996l) . 

The advent of space telescopes, combined with ground- 
based observatories, has allowed astronomers to observe vapours 
and ices of water in several phases of star formation. Water 
vapour is present in cold molecular clouds and prestellar cores 
but only with low abun dances (X gas (H20) ~ 10~ 8 - 10~ 9 rel- 
ative to H nuclei, see Bergin & Snell 120021: iKlotz et all l2008t 
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ICaselli et"aT1 120 1 Oh . In cold clouds, w ater is likely condensed 
in ices (X lV?g (H 7 ,0) ~ 5 x 10~ 5 - 10~ 4 . IWhittet & Puled fl99ll: 
iPontoppidan et al.l [2004>. Hot corinos and outflows of low-mass 
Class protostars show higher abundances of g as phase water, 
with a bundances of a few 10~ 6 in hot corin os ( Ceccarelli et all 
2000; ICoutens et ail 120121: iKristensen etafll2012h . whilst pro- 
tostar outflows show higher abundan c es up to a few 10~ 5 
dLiseau et alj|1996t ILefloch et all [20101: IKristensen et alJfeOlQl 
2012). Water ice has also been observed towards cold proto- 
stellar en velopes with similar abundances to mo lecular clouds 
(~ 10~ 4 . IPontoppidan et al.ll20Q4T: iBoogert et al.ll2008l) . Recent 
infrared observations have also shown the presenc e of water in 
proto planetary disks in different states: w ater ice dTerada et all 
1 20071) and hot and co ld water vapour (ICarr & Najital 12008: 
Hogerheii de et al.l201 ll with abundances of 10 -4 and lower than 
10~ 7 , respectively). Analysis of debris di sks show that dust par- 
ticles are covered by ice mixtures (e.g. iLi & Greenbergl ll998: 
iLebreton et aDl2012h . 

To summarise, it is now accepted that water is present during 
all phases of the star formation process. However, its evolution 
from molecular cloud to planetary system still remains unclear. 
The deuterium fractionation can help us constrain its formation 
and its evolution. First, it can probe the formation pathways of 
water observed in the early stages because of its sensitivity to 
the physical conditions. Second, it allows us to investigate its re- 
processing in protoplanetary disks, and eventually to determine 
whether water on Earth has an interstellar origin. Comparing the 
HDO/H2O ratio in comets and Earth is, for example, important 
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for evaluating the possible co ntribution of comets for transfer- 
ring water in Earth's oceans (lOwen & Bar-Nunlll995l) . Recent 
Herschel observations have reported a D/H ratio of water (0.014 
%) in the Jupiter family comet 103P/Hartley2 originating in the 
Kuiper Belt, very similar to the value for the Earth's oceans, sup- 
porting the hypothes is that a part of water comes from comets 
(iHartoghet allhoill) . 

The HDO/H2O ratio has recently been evaluated in 
the gas phase of low-mass Class protostar envelopes 
with val ues varying from less than 1 0~ 4 in NGC1333- 
IRAS4B (J0rgensen_& van Dishoeck 2010) to more than 10 
in NGC1333-IRAS2A (|Uu_eiaLj|2011|)- The low-mass Class 
protostar IRAS 16293-2422 seems to have the most reliable 
value since the main isotopologue (via H^O) and its simply 
and doubly deuterated isotopologues have been observed several 
times via ground-based and space telescopes (IC eccarelli etall 
2000: lPariseetal.ll2005t iButner et all 120071: IVastel et al.ll2010t 
Coutens et al. 2012). The most recent work by Cout ens et al.l 
(1201 2l) reports an HDO/H 2 ratio of ~ 3 % in the hot corino, 
~ 0.5 % in the cold envelope, and ~ 5 % in the photodes- 
orption layer of the for eground cloud. Obs ervations of D2O by 
IButner et al.l (120071) and lVastel etaPfcOlOl) give a D 2 0/H 2 ra- 
tio of ~ 10~ 3 . Using the 4.1 OD stretch band, solid HDO has 
been observed towa rds a sample of low- an d high-mass stars by 
(iDartois et al.l l2003) and Paris e et aTI(l2003l) . These observations 
provide an upper limit for the solid abundance ratio between 0.2 
and 2 %. Although these values are much highe r than the cos - 
mic elemental abundance of deuterium (1.5x 10~ 5 lLinskvl2 003). 
water seems to be less deuterated than the other molecules also 
mainly formed on grain surfaces, such as formaldehyde and 
methanol. Figure Q] graphically shows that water possesses a 
lower level of deuteration, with a mean HDO/H2O ratio of ~ 3 
%, whereas other molecules have a fractionation of their singly 
deuterated isotopologue higher than 10 %. 

H 2 0, but also H 2 S, H 2 CO and CH 3 OH (and HCOOCH3) are 
thought to be mainly formed on grain surfaces via the hydro- 
genation of accreted O (or 2 ), S, and CO from the gas phase. 
Cold gas-phase che mistry produces wat er vapour with low abun- 
dances (a few 10~ 7 : iBergin et alJ uOOO). whilst endothermic re- 
actions producing warm water are efficient at T > 250 K. In 
contrast, water ice desorbs into the gas phase in the envelope of 
Class protostars in large quantities at T ~ 100 K. Therefore, 
the deuteration of water measured by millimetric observations 
likely reflects the deuteration of its icy precursor governed by 
the accretion of gas phase H and D atoms. Indeed, the timescale 
needed to significantly alter the deuteration after evaporation 
in warm gas is longer than the typical age of Class proto- 
stars ( ~ 10 5 versus ~ 10 4 y r, Icharnlev et al.lll997b lAndre et al.l 
2000). iRoberts et all (120031) showed that the gas phase atomic 
D/H ratio increases with the CO depletion. Molecules that are 
formed in the earlier stages of star formation, when the CO de- 
pletion is low, would, therefore, show lowe r deuteration. Based 
on our astrochemical model GRAINOBLE ( Taquet et al.ll2012al 
hereafter TCK12a), we theoretically confirmed this hypothe- 
sis by successfully repr oducing the observe d formaldehyde and 
methanol deuterations dTaquet et al.ll2012bl hereafter TCK12b). 
The difference in deuteration between these two molecules is 
explained by the earlier formation of formaldehyde compared 
to methanol, when the D/H is lower. We also demonstrated 
the necessity t o introduce the abstraction r eactions exp e rimen - 
tally sho wn bvlNagaoka et all (120071) and iHidaka et alJ (l2009h . 
Similarly, ICazaux et al.l (1201 ll) showed the possibility that wa- 
ter is formed through reactions involving H 2 . According to these 
authors, instead of reflecting the atomic D/H ratio, water deuter- 
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Fig. 1: Deuterium fractionation of several molecules assumed to 
be partly (or mainly) formed on interstellar grain surfaces and 
observed around low-mass Class protostars. Blue triangles: 
simple deuteration; red diamonds: double deuteration; green 
squares: triple deuteration. 

References. 

(a) Cosmic D/H: lLinskyl d2003l) 
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W H2Q : HDO/H 2 bv lCoutens efaU J20T21) . D 2 0/H 2 bv lVastel et all 
fcOlOh towards IRAS 16293. 

<c) H2S: HD S/H,S bvlvan Dishoecketai] fl99l) towards IRAS 16293, 

D 2 S/H 2 S bv lVastel et al] d2003h towards I RAS4A. 

w NH: ND/NH bv lBacmann et al] <201Ch towa r ds IRA S 16293. 

(e) NH3: NH 2 D/NH7 ~bv Ivan Dishoeck et al] d 19951) towards IRAS 

16293; NHD ,/NH, bvlLoinard et al] d2001l) towards IRAS 16293, and 

ND1/NH3 bv lvan der Tak et al] d2002|) towards IRAS 4A . 

W H2CO: H DCO/H 2 CO by Ivan Dishoeck et al] JT9951) towards IRAS 

16293 and by lParise _etalJ 0006]) toward s sev en low-mass pro t ostars ; 

D 2 CO/H 2 CO by ICeccarelli et al] d 19981) and ICeccarelli et al] d2001h 

toward s IRAS 16293 and bv lParise et aU d2006h and lRoberts & Miiiail 

(2007) towards a sample of low-mass protostars. 

& CH 3QH: CH 2 DOH/CH 3 OH and CHD 2 OH/CH 3 OH bv lParise etall 
( 2006) towards seven low-mass protostars, CD3OH/CH3OH by 
Parise et al] d2004 towards IRAS 16293. 

m HCOOCH3: DCOOCH 3 /HCOOCH 3 by iDemvk et al] J20lch to- 
wards IRAS 16293. 



ation should scale with the gas phase HD/H 2 ratio (~ 10~ 5 ) at 
low temperatures. 

The linear relationship between the observed column den- 
sity of water ice and the vi sual extinction , abov e the threshold 
of A v ~ 3.2 mag found by Whi ttet et all (Il988l) . suggests that 
water ice starts to form significantly at low visual extinctions. 
Other solid compounds are believed to form along with water 
ice. The linear relationship between carbon dioxide and water 
column densities, with a column density ratio Af(C0 2 )//V(H 2 0) 
of 18 %, suggests that these two molecules form in parallel 
(Whittet et al. 2007). This conclusion is supported by compar- 
isons between observed and laboratory band profiles showing 
that C0 2 is mainly located in a polar water-rich mixture whilst 
a non-polar component (pure C0 2 , or CO:C0 2 ) exist in very 
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low quantities dGerakines et al.lll999b iPontoppidan et al.l l2008). 
Other solid organic compounds, such as CO or CH3OH, have 
also been observed in large quantities (~ 30 % for CO, and up 
to 20% for methanol) but at higher visual e xtinctions, above Ay 
thresh olds of 8-9 and 15 mag respectively (IWhittet et al.ll2007L 
1201 lb . Consequently, one needs to include the formation of all 
other solid molecules in order to correctly study the formation 
of deuterated water ice. 

In this article, we extend our study to the formation 
and the deuteration of the interstellar water ice formed in 
molecular clouds and in prestellar cores, using the multilayer 
GRAINOBLE model. Our goal is to explore the influence of 
the physical conditions (representative of different typical cloud 
stages) on ice formation and water deuteration. This is the first 
time that such a systematic study has been done. Besides that, 
it includes the crucial influence of the H2 ortho/para ratio on the 
water deuteration. The astrochemical model is presented in Sect. 
12 In Sect. [3] we study the formation of typical grain mantles 
and show the influence of several physical and chemical param- 
eters on the deuteration of water ice. In each section, we sum- 
marise the main ideas with headings and concluding remarks. In 
Sect. |U we compare our with previous model predictions and in 
Sect. [5] with published observations of water, formaldehyde, and 
methanol. 

2. Modelling 

2.1. Overview of the GRAINOBLE model 

GRAINOBLE is a gas-grain astroche mical model based on th e 
rate equations approach introduced by Hasega wa et al.l dl 992b . 
A detailed presentation of GRAINOBLE can be found in 
TCK12a. Briefly, GRAINOBLE couples gas-phase and grain- 
surface chemistry. In total, our chemical network consists of 341 
(gaseous and solid) species and 3860 (gas phase and grain sur- 
face) reactions. The gas phase chemistry is described in detail in 
section 12.21 The grain surface chemistry processes are the fol- 
lowing: 

i) The accretion of gas phase species onto the grain surfaces, as- 
sumed to be spherical. 

ii) The diffusion of adsorbed species via thermal hopping. 

iii) The reaction between two particles via the Langmuir- 
Hinshelwood mechanism, once they meet in the same site. The 
reaction rate is the product of the number of times that the two 
reactants meet each other and the transmission probability P r of 
reaction. 

iv) The desorption of adsorbed species into the gas phase via 
several processes: 

- thermal desorption; 

- cosmic -ray induced heating of grains; 

- chemical desorption caused by the energy release of exother- 
mic reactions. 

This last process is an up grade with respect to TCK12a. 
Following iGarrod et al. (2007p, we assume a value of 0.012 for 
the factor a (the ratio of the surface-molecule bond frequency to 
the frequency at which energy is lost to the grain surface) since 
it seems to be the m ost consistent value given by molecular dy- 
namics simulations (Kroes & Andersson 2005). 
In addition, in the present work, we take the effect of the UV 
photolysis on the ices into account (see section 12.5b . because 
this is important for the formation of H2O at low visual ex- 
tinctions. As suggested by l aboratory experiments, cold man- 
tle bulks are mostly inert ( see lWata nabe et al. 20 041 iFuchs et al.l 
l2009t lloppolo et al1l2010l) . Therefore, we follow the formation 



of grain mantles with a multilayer approach in which the outer- 
most only layer is reactive, whilst the mantle bulk remains inert 
(see TCK12a for more details). 

2.2. Gas-phase chemical network 

We consi der the gas-phase che mical network from the KIDA 
database ( Wakelam et al.l 120 1 2h for seven elements: H, He, C, 
N, O, S, and Fe, giving a total of 258 gaseous species. Sulphur 
and iron are introduced to consistently study ion chemistry. 
Reactions involving atomic Fe and S play a significant role in 
the destruction of H + and H3 whilst S + , and to a lesser extent 
Fe + , is believed to be one of the mo st abundant ions before the 
CO depletion (see lFlower et al.ll2005l) . Relative to the KIDA net- 
work, we modify the rate coefficient o f the cosmic-ray disso - 
ciation of H2 yielding H+H. Following Dalgar no et alj (fl999). 
we assume y - 0.5, where the rate kn ss (in s _1 ) of cosmic ray 
dissociation reactions is given by kdiss - 7% (<T is me cosmic- 
ray ionization rate). Before ice formation, bare grains are con- 
sidered. The recombination efficiency of atomic hydrogen and 
deuteriu m are assumed to be unit y, fo llowing the theoretica l 
works by ICazaux & Tielensl (|2004|) and ICuppen et ail (l2010bl) . 
who considered chemisorption interactions. The initial elemen- 
tal abundances in the gas phase consi dered in this wo r k are listed 
in Table [Hand follows the work by Wakel am et al.1 (1201 Ol) and 
Linskv (2003). To model the water formation at low visual ex- 
tinctions (see Introduction), we consider the depth-dependent 
self-shielding of H2, HP, a nd CO, using the Meudon PDR code 
dLe Petit et alJl2002ll2006l) . 

Since standard gas-phase models overpredi ct the O2 abun- 
dance with respect to w hat is observed (see iGoldsmith et al.l 
1201 lb iLiseau et al.ll2012b we also ran a grid of models by arti- 
ficially decreasing the formation rate of O2 by a factor of 10. 
We checked that the [02]/[0] abundance ratio remains indeed 
ten times lower than the "standard" case throughout the whole 
calculation. 

Deuteration of water and other molecules formed on inter- 
stellar grains strongly depend on the gas phase abundances of 
H, D, H2, HD, and D2. In turn, the abundance of these gaseous 
species, and more particularly D, depends mainly on the deuter- 
ation of Hi which is function of the degree of CO and Ng 
freeze-outs onto the grain surfaces (see Roberts & Millar! 120001 
Rober ts et al 120031 The deuterium gas phase chemical network 
is based on that of TCK12b. The major differences with respect 
to TCK12b are the inclusion of the ortho-to-para (hereinafter 
opr) H2 ratio and reactions involving N2. 

Reactions between Hi isotopologues and H2, reducing the 
deuterium fractionation, are endothermic. Therefore, these reac- 
tions cannot occur in cold-cloud conditions if para-H2 only is 
considered. However, H2 is also believed to exist in the ortho 
spin state, higher in energy (170 K), since H2 probably forms 
on grains with an ortho-to-para ratio (opr) of 3. Reactions be- 
tween ortho-H2 and ortho-H2D + (ortho-HDJ, ortho-Dp can re- 
duce the degree of deu t eration of H^ significantly at low tem- 
peratures. Flow er et al.l ([2006) have shown that the fractiona- 
tion of H3 , hence, the abundance of atomic D are strongly re- 
duced when the H2 opr is higher than 10~ 4 . Consequently, we 
enlarge the chemical network, relative to TCK12b, by consid- 
ering the H3-H2 system whose new reactio n rate coefficients 
have been computed by Hugo et al. (2009). Ion-neutral reac- 
tions between H3 isotopologues (including their different spin 
states) and CO or N2, electronic recombinations, and recombina- 
tions on electronegative charged grains have also been included 
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following iRoberts & MiUaM 2000): Robe rts et al .1 (120031 12004): 
IWalmslev et alJ (120041) and Pagani et all (120091) . 

The actual opr of H2 in molecular clouds is still highly un- 
certain. The initial value of H2 opr formed on grain surfaces 
is most likely 3, as re cently confirmed by the experiment by 
IWatanabe et al] d2010h conducted on amorphous solid water. 
Proton-exchange reactions in the gas phase would then con- 
vert ortho-H2 to para-H2, decreasing the o pr of H? towards th e 
Boltzmann value (~ 3 x 10~ 7 at 10 K, see iFlower el al. 2006). 
Recent experimental studies have also demonstrated that H2 un- 
dergoes a nuclear spin conversion from th e ortho to the para spin 
state on amorphous solid w a ter (ASW) (ISugimoto & Fukutanil 
l201lt IChehrouri et al.ll20TTt lHama et al.ll2012l) . The influence 
of the H2 o pr on absorption lines o f formaldehyde has been 
observed by iTroscompt et al.1 (120091) who deduce that the H2 
opr is much lower than 1. Indirect estimates of H2 based on 
the compariso n with chemical models suggest valu es of about 
10 3 - 10~ 2 dPagani et alj|2009t iDislaire et al]|2012l) . Given the 
relative uncertainty in this value and its importance in the molec- 
ular deuteration process, in this work we assume the H2 opr as a 
free parameter constant with time. 



Table 1: Initial elemental abundan ces in the ga s phase with re- 
spect to hydrogen nuclei (from iLinskvl 12003: Wakelam et al.l 
2010). 



Species 


Abundance 


H 2 


0.5 


HD 


1.6 x 1(T 5 


He 


0.09 


C 


1.2 x 10~ 4 


N 


7.6 x 10~ 5 


O 


2.6 x 10~ 4 


S 


8.0 x 10~ 8 


Fe 


1.5 x 10~ 8 



2.3. Chemical network on grain surfaces 

2.3.1 . Formation and deuteration of water ice 

We consider a chemi cal network based on the work by 
iTielens & Hagenl d!982l) modified following the results of sev- 
eral recent experimental works as described below. The simplest 
formation pathway towards solid water is the sequential hydro- 
genation of atomic oxygen: 



O + H — > OH 



OH + H — > H 2 0. 



(1) 

(2) 



This reaction channel was experimentally measured to oc- 
cur in cold condi t ions, probably via barrierless re actions 
dHiraoka et aUll998HDuiieu et al]l201rtljing et al.ll201 ll) . 

In addition, water can also be formed from different channels 
involving O2 or O3 as follows. First, water can be formed from 
molecular oxygen, either from the gas phase or formed on grains, 
following the reaction channels: 



O2 + H — > H0 2 
H0 2 + H -> H 2 2 
H 2 2 + H -> H z O + OH. 



(3) 



These pathways have been experimentally de monstrated by 
Miva uchi etafl (120081) : lloppolo et al.l (l2010l) and ICuppen et al.l 
d2010al) at temperatures of 10 K. iMivauchi et alJ d2008l) also 
showed an isotope effect in the formation of water from hydro- 
gen peroxide, implying the possibility of tunnelling through an 
activation barrier for this reaction. 

Sec ond, water can be formed from ozone on interstellar 
grains (Cuppen & Herbst 20Q2t iTaquet et al.l 120 1 2ah following 
the reaction 



O3 + H -> 2 + OH 



(6) 



Then, O2 and OH can continue to react to form w ater 
via the reactions described above. Mokrane et al .] (120091) and 
iRomanzin et al.l d201 ll) experimentally showed the efficiency of 
this reaction by observing the presence of water (HDO, D2O) 
after the irradiation of solid O3 on water ice by H (D) atoms. We 
treat this reaction as barrierless. 

With their micro scopic Monte Carlo model, 
Cupp en & Herbstl d2007l) have concluded that molecular 
hydrogen plays a key role in the formation of water ice in 
molecular clouds. When hydrogen is mostly in its molecular 
form, water is formed at ~ 70% by the reaction 



OH + H 2 -> H 2 + H. 



(7) 



lOba et all d2012l) experimentally showed that HDO formation 
from OH is ten times less efficient than the formation of H2O 
from OH, implying an isotope effect and therefore the possibil- 
ity of tunnelling through an activation barrier in this reaction. 

OH radicals can also recombine if the ice temperature is 
high enough (4 K in their experiment) to allow their mobility. 
lOba et all (120 1 ll) experimentally determined branching ratios for 
the reactions 

OH + OH -> H z O + O (8) 

OH + OH -> H 2 2 (9) 

of 0.2 and 0.8, respective ly. 

Unlike ICazaux et al] (120101 1201 lb , we do not include the 
O+H2 reaction since it is unlikely to proceed at low temperatures 
(10-20 K) given its hig h endothermicity (960 K; iBaulch et al] 
1992). IOba et al]d2012l) confirmed that the co-deposition of cold 
O atoms with H2 at 10 K does not result in the formation of water 
but only of O2. 

We enlarge the water-formation network by including the 
deuterated counterparts of all reactions mentioned above. Due to 
the higher mass of deuterated species with respect to their main 
isotopologue, the reaction rates of barrierless reactions involv- 
ing deuterated species are decreased. A careful treatment of the 
transmission probabilities (i.e. probability of tunnelling through 
the activation barrier) of all reactions possessing a barrier is de- 
scribed in section l2~4l 



2.3.2. Formation and deuteration of other ices 

The accretion of CO and O particles onto interstellar grains can 
lead to the formatio n of carbon dioxide (CO2) mainly via three 
reaction channels (Ruffle & Herbst 2001) 



HCO + O 



(4) CO+O- 

(5) CO + OH 



-» CO2 + H 
C0 2 

-> CO2 + H 



(10) 

(ID 
(12) 
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Gas phase experiments have shown that reaction (fTOt is barrier- 
less (e.g. Bau lch et alJ l2005) although it has never been studied 
on interstellar ice analogues so far. 

Reactio n (fTTt is thought to have a high activation energy 
dTalbi et al.l 120061) . Laboratory experiments have shown that 
solid carbon dioxide can be formed from this reaction with low 
efficiency, at least two or ders of magnitude lower than astro nom- 
ical observations dRoser et al.1 12001; Raut & Baragiola 201 ll) . 

Formation of solid carbon dioxide, from OH radicals and CO 
molec ules, has been observed even at very low temperatures (10 - 
20 K, lOba et alJl2010t lloppolo et alJl201 it iNoble etalJl201lh . 
However, the exact path way of the formation of CO2 is still un- 
certain. Following lOba et all (1201 Ol) who observed a weak band 
attributed to the HOCO radical and other theoretical works, we 
propose a reaction pathway for the CO2 formation from CO and 
OH, as described in Appendix A.2. 

The chemical network presented in TCK12b, with the same 
relative rates, is used to study the formation and the deuteration 
of formaldehyde and methanol. We also consider the formation 
of deuterated methane and ammonia from barrierless addition 
reactions of solid atomic carbon and nitrogen. A full list of sur- 
face reactions considered in this work is presented in Appendix 
B. 



2.4. Eckart model and reaction probabilities 

In previous gas-grain astrochemical models, the transmission 
probability of exothermic surface reactions has been approxi- 
mated by the exponential portion of the quantum mechanical 
probability for tunnelling through a square potential barrier of 
width a. However, this approach has two main limitations: 

i) it does not fit the potential energy profile of reactions correctly; 

ii) the width a is unknown, although in most astrochemical 
models it is arbitrari l y fixed to 1A (iTielens & Hagenl Il982l 
lHasegawaetalJl!992l) . iGarrod & Paulvl (1201 ll) reproduced the 
value of the transmission probability of a few reactions with 
square potential barriers and deduced a width of about 2 A. 
However, the deduction of this width is based on poorly con- 
strained activation barriers. 

Since square potential barriers do not allow us to accurately 
estimate the values of the transmission probabilities, we com- 
pute the trans mission probabilities of all the reactions u sing the 
Eckart model (Eckart| [l930t iJohnston & Heicklenll 19621) . which 
fits an approximate potential energy surface. A full description 
of the Eckart model, our quantum chemistry calculations and 
the computations of the transmission probabilities are given in 
Appendix A. 

Briefly, the formation of deuterated water includes two reac- 
tion channels which have an activation barrier: reactions dH) and 
©. Reactio n (|5]l has been theoreti cally studied by Koussa et al. 
(2006) and Ellings on et al.l ([2007) but only for the main iso- 



depends on the area of the potential energy profile, the Eckart 
model gives higher probability of transmission than the square 
barrier (1.4 x 10~ 7 versus 1.2 x 10~ 8 ). 



topologues. Therefore, quantum chemistry calculations are con- 
ducted for all the reactions (|5} including deuterated isotopo- 
logues (see Appendix A for more details on the calculations). 
Concerning reaction (|7}, the transmission probabilities of all the 
deuterat ed counterparts have been deduced from the theoretical 
work bv lNguven et al.l d201 ll) . 

Figure 12 shows the Eckart and the symmetric square poten- 
tials as a function of reaction coordinates for the reaction <(3j 
computed in this work. As suggested previously, it can be seen 
that the Eckart potential is far to be symmetric. Furthermore, the 
reaction profile is thinner than the square barrier. Since the quan- 
tum tunnelling probability of transmission through this barrier 




-1 1 

Reaction coordinate x 

Fig. 2: Potential energy profile as function of the reaction coor- 
dinate of the H2O2 + H — > H2O + OH reaction computed from 
ab-initio calculations (black solid curve) and adopting a sym- 
metric square barrier of a width of 1 A (red dashed curve). See 
Appendix A for more details on calculations. 

Table [2] lists the reactions with an activation barrier, as well 
as the input parameters needed for computing the transmission 
probabilities. It also compares the transmission probability com- 
puted with the Eckart model and with a symmetric square poten- 
tial barrier of the same activation energy and adopting a width 
of 1 A (the value commonly used in most gas-grain astrochemi- 
cal models). The comparisons between the two approaches show 
that the assumption of a square barrier width of 1 A tends to un- 
derestimate the reaction probabilities for all the reactions by up 
to seven orders of magnitude. 

The transmission probability of the CO+H reaction com- 
puted with the Eckart model is in good agreement with the range 
of values we deduced in TCK12a. In this latter work, we varied 
the transmission probability P r of this reaction and found that a 
transmission probability higher than 2 x 10~ 7 was needed to re- 
produce the solid CH3OH/CO ratio observed towards high-mass 
protostars. 

2.5. Photodissociation and photodesorption of ices 

In addition to Langmuir-Hinshelwood chemical reactions, 
we also consider the effect of FUV (6 - 13.6 eV) pho- 
tons on ices, following the res ults of molecula r dyn amics 
(MP) simulations carried out by lAndersson et al.l (120061) and 
[And ersson & van Dishoeckl (120081) . They showed that an amor- 
phous ice absorbs UV photons in a 1-2 eV narrow band peaked 
at ~ 8.5 eV, in g ood agreement with experimental works by 
Kobavashij (1 19831) . reaching a maximal absorption probability 
of 7 x 10 3 per monolayer. We convolve the absorption spec- 
trum with the emission spectrum of the in terstellar radiation 
field (ISRF) deduced by Mat his et alJ (l983) and with the emis- 
sion spectrum of H2 excited by the secondary electrons pro- 
duced by the cosmic - rays i onization of H2 (CRH2RF) com- 
puted bv iGredel et al.l (0.987). The absorption probabilities P a t s 
of each monolayer integrated along the 6 - 13.6 eV band are 
therefore equal to 1.51 x 10 3 for ISRF and 1.03 x 10 3 for 
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Table 2: List of reactions having a barrier, with the forward, backward activation barriers, the imaginary frequency of the transition 
state, the transmission probabilities computed with the Eckart model and the square barrier method (using a barrier width of 1 A). 



Reactions 














Probability type 


V f (K) 


Vj (A) 


Vs 


'r.Eckan 


-* r,$quare 


Reference 


OH 


+ 


H 2 


"~ * 


H 2 


+ 


H 


Eckart 


2935 


10209 


-1293 


4.07(-07) 


1 .49(- 13) 


1 


OH 


+ 


HD 


~~ * 


H 2 


+ 


D 


Eckart 


2855 


9396 


-1259 


3.62(-07) 


7.91(-16) 


1 


OH 


+ 


HD 


"~* 


HDO 


+ 


H 


Eckart 


3051 


10508 


-970.4 


1.00(-09) 


2.44(-16) 


1 


OH 


+ 


D 2 


~* 


HDO 


+ 


D 


Eckart 


3026 


9556 


-955.3 


8.07(-10) 


8.02(-18) 


1 


OD 


+ 


Hi 


~* 


HDO 


+ 


H 


Eckart 


2789 


10246 


-1293 


8.74(-07) 


2.88(-13) 


1 


OD 


+ 


HD 




D 2 


+ 


H 


Eckart 


2900 


10871 


-970.1 


2.81(-09) 


5.2(-16) 


1 


OD 


+ 


HD 


— * 


HDO 


+ 


D 


Eckart 


2703 


9736 


-1258 


7.99(-07) 


1.76(-15) 


1 


<JL> 


+ 


D 2 


— * 


D 2 U 


+ 


u 


Eckart 


Z5 /U 




1 

-yjj.L 


OA/" 






H 2 2 


+ 


H 


— * 


H 2 


+ 


OH 


Eckart 


2508 


36358 


-1054 


1.37(-07) 


L18(-08) 


2 


1 1 (~\ 


+ 


D 




tlDV 


+ 


1 it 1 

un 


Eckart 


Zdjj 


V1 1 1 Q 
J 1 1 to 


-o4Ji. / 




o.oj(-1Z; 


z 


I I FA / A 

rlDU? 


+ 


1 1 

H 




HDU 


+ 


UH 


Eckart 




ibliy 


1 n?i 


1 IK f\H\ 


1 n/ no\ 
1.1 /(-UoJ 


Z 




+ 


1 1 
rl 




H 2 U 


+ 


UL) 


Eckart 


zjZ4 


1 AAAQ 


-1U3 i 


1 OO/ fY7\ 
L.ZZ(-\J 1 ) 


1.1 /{-\JO) 


Z 


riDLb 


+ 


JJ 




HLHJ 


+ 


UL) 


Eckart 


liby 


jboZZ 


O /I A A 

-o4o.O 


j.Zo{-\jy) 


Q AA/ 1 T\ 


Z 


I I FA / \ 

riJJLh 


+ 


JJ 




D 2 U 


+ 


Url 


Eckart 


lib! 




Q A A 1 

-840.1 




O.OD(-1ZJ 


Z 


U2U2 


+ 


rl 




HDU 


+ 


UL) 


Eckart 


234U 




-1U32 


1 HO/ A r 7\ 

l.Uo(-U / ) 


1.1 /(-Uo) 


Z 


D2U2 


+ 


D 




D 2 U 


+ 


UL) 


Eckart 




ib Hi 


O/IO 

-B4z.y 


a to/ nn\ 


/i n/ 1 1\ 


z 


CO 


+ 


H 


— * 


HCO 






Eckart 


1979 


8910 


-793.6 


1 no/ rn\ 

1.92(-07) 


1 O "} / flO\ 

1.83(-08) 


3 


CO 


+ 


D 




DCO 






Experiments 








1.92(-0o) 


1.83(-09) 


4 


H 2 CO 


+ 


H 




CH3O 






Experiments 








9.60(-08) 


l~\ ~\ C t f\l \ \ 

9.15(-09) 


4 


H2CO 


+ 


D 




CH 2 DO 






Experiments 








r\ /~ f\/ r\c\\ 

9.60(-09) 


9.15(-10) 


5 


H 2 CO 


+ 


D 


— * 


HCO 


+ 


HD 


Experiments 








9.31(-08) 


00/ r\f\\ 

8.88(-09) 


4 


H2CU 


+ 


L) 


— * 


HULL) 


+ 





Experiments 








9.31(-Uoj 


00/ new 


4 




+ 


O 

H 


— * 


LH 2 DO 






Experiments 








1 11/" m\ 
l.ll(-0/) 


1 r\£.{ r\c\\ 

i.Uo(-i)y) 


4 


HDCO 


+ 


D 


— * 


CHD 2 






bxpenments 








1.1 1 (-07) 


1 r\£~{ r\c\\ 


5 


HDCL) 


+ 









+ 


irtU 


Experiments 








1.34(-U/) 


1.4o(-Uo) 


c 

j 


HDCO 


+ 


D 




DCO 


+ 


HD 


bxpenments 








(i ti / r\o\ 


o.87(-U9) 


5 


HDLU 


+ 


L) 






+ 





Experiments 








n 11/ no\ 

9.31(-Uo) 


0.0 /(-UV; 


c 

J 




+ 


u 
1 1 




pt4,r»o 






Experiments 








1 97/' 07^ 


1 7 1 /■ ns^ 

1 .Zl^-UO,) 


A 


D2C0 


+ 


D 




CD30 






Experiments 








1.27(-08) 


1.21(-09) 


5 


d 2 co 


+ 


H 




DCO 


+ 


HD 


Experiments 








7.30(-08) 


6.95(-09) 


4 


CH,OH 


+ 


D 




CH 2 OH 


+ 


HD 


Experiments 








2.88(-07) 


2.75(-08) 


4 


CH 2 DOH 


+ 


D 




CHDOH 


+ 


HD 


Experiments 








1.92(-07) 


1.83(-08) 


4 


CHD 2 OH 


+ 


D 




CD 2 OH 


+ 


HD 


Experiments 








1.50(-07) 


1.43(-08) 


4 


CH,OD 


+ 


D 




CH 2 OD 


+ 


HD 


Experiments 








2.88(-07) 


2.75(-08) 


5 


CHoDOD 


+ 


D 




CHDOD 


+ 


HD 


Experiments 








1.92(-07) 


1.83(-08) 


5 


CHD 2 OD 


+ 


D 




CD 2 OD 


+ 


HD 


Experiments 








1.50(-07) 


1.43(-08) 


5 



Notes. The probability type refers to the method used for computing the probability. "Experimen ts" means that the transmission probability 
is deduced from the CO+H reaction by considering the relative rates experimentally measured by Hida ka et al.l d2009h (see TCK12b for more 
details). 

References. l: lNguven et"aTJ ( l2011l) . 2: This work. 3: lPeters et a0(l2012h . 4: lHidakaet all<2009l) . 5: lTaauet et aTU2012bl) . 



CRH2RF. lAndersson & van" Dishoeckl (120081) also showed that 
the photofragments, OH and H, display different trajectories 
(desorption, trapping, or mobility on the surface) and they com- 
puted the yield of each trajectory as fu nction of the mo n olayer 
Subse quen t experimen t al studies b y Yabush ita et al.l (2006, 
120091) and lHama et alJ (120091 l2Q10h have confirmed the dif- 
ferent photofragment trajectories revealed by MD simulations. 
However, our model only considers one chemically reactive 
layer (the outermost). We deduce th e yield of each trajectory by 
averag ing the yields computed by lAndersson & van Dishoeckl 
(2008) on each layer. The list of trajectories as well as their aver- 
aged yield Y is given in Table [3] The rate of each trajectory can 
be deduced, via the following equation 

RphotoAs' 1 ) = Yi x F uv x cr(a d ) x 5 x P abs x l/N s (13) 

where Fuv (cm -2 s" 1 ) is the photon flux of the 6 - 13.6 eV UV 
band which is absorbed by the amorphous ice, cr(a d ) (cm 2 ) is 
the cross section of interstellar grains (= 7r(a<//2) 2 where a,/ is 
the grain diameter), 5 refers to the absorption of the five outer- 
most layers, and Af, is the number of sites on the grain surface. 



Owing to the lack of quantitative data on the photodesorption of 
atoms and on the photodissociation of hydrogenated (and deuter- 
ated) molecules other than water on ASW ice, we consider the 
same absorption probability for all atoms and molecules as for 
water. For atoms, the desorption probability upon absorption is 
assumed to be unity, whilst we consider the same outcome prob- 
abilities for hydrogenated molecules. 



Table 3: List of trajectories after H2O photodissociation with 
their probability (deduced from Anders son & van Dishoeckl 
120081 see text). 



Outcome Probability 



H 2 0, t , 




Hgas + OHfc e 


0.5 


H 2 0, ce 




Hfcg 4- OH, C(J 


0.2 


H 2 0,« 




Hfc e + OH gfl! 


2.2 x 10~ 3 


H 2 i ce 




Hgas "t" OHj3 ( y ? 


6.8 x 10~ 3 


H 2 0,„ 




H 2 O gfl! 


4.0 x 10~ 3 


H 2 0;„ 




H 2 0, <:( , 


0.28 
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We follow the experimental results by iFavolle et al.l (1201 ll) 
for wavelength-dependent CO photodesorption between 7.5 and 
13 eV. The convolution of the CO photodesorption spectrum 
with the ISRF and CRH2RF fields gives integrated photodes- 
orption yields of 1.2 x 10~ 2 and 9.4 x 10~ 3 photon 1 molecule 
_1 , respectively. Owing to the lack of data on the wavelength- 
dependent photodesorption of other molecules, we consider the 
same photodesorption rates for O2, and N2. 

2.6. Binding energies 

Comparisons between the observed absorption 3 fim band of wa- 
ter and laboratory experiments have shown that grain mantles are 
mainly composed of high-density amorpho us solid water (ASW) 
dSmifh et alj|1989t Uenniskens et aliri995l) . Therefore, binding 
energies of adsorbed species relative to ASW must be consid- 
ered. 

It is now accepted that light particles (H, D, H2, HD, 
D2) show a distribution of their binding energies relative to 
ASW, depending on the ice properties, the adsorpt ion condi- 
tions, and the coverage of the accreted particles. iPerets etail 
have experimentally highlighted the influence of the ice 
density on the HD and D2 binding energies. Molecules ad- 
sorbed on low-density amorphous ices (LDI) desorb follow- 
ing three desorption peaks that are at lower temperatures than 
the single observed desorption peak of molecules evaporating 
from a high-density amorpho us ice (HDI). By depositing D2 on 
ASW. lHornekasr et all d2005) showed that D 2 is more efficiently 
bound to porous surfaces. Furthermore, binding energies follow 
broad distributions between 300 and 500 K and between 400 
and 600 K fo r non-p orous and porous ASW ices, respectively. 
Ami aud et al.1 d2006l) studied the link between the binding en- 
ergy distribution of D2 with its coverage on a porous ASW ice 
and with ice temperature. More particularly, they showed that 
the D2 binding energy decreases with H2 coverage, from ~ 700 
K to ~ 350 K whilst distributions broaden. 

The experimental results have been supported by theo reti- 
cal MP calculations. For exampl e, iHornekaer et al.l (120051) and 
lAl-H alabi & v an Dishoeckl d2007l) found that the binding energy 
distributions are essentially a consequence of the variation in 
the number of water molecules surrounding the adsorbed par- 
ticle. Thus, the binding energy of H shows a broader distribution 
peaked at higher values on porous and irregular amorphous wa- 
ter ice than on a structured crystalline ice. 

It is believed that deuterated species are more efficiently 
bound with i ces than their main isotopologues because of their 
higher mass (lTielensl [l983). However, the difference in binding 
energy between H, H2, and their deuterated counterparts s till re- 
mains poorly constrained. Experiments by IPerets et aD d2005h 
and iKristensen et al.l d2011l) have shown that distribution peaks 
of the H2, HD, and D2 binding energies are very close (less than 
5 meV ~ 60 K). These differences remain much smaller than 
the typical full-width-at-half-maximum of binding energy distri- 
butio ns shown by MP simulations (~ 195 K for a amorphous 
ice lAl-Halabi & van Dishoeckl l2007). we, therefore, assume the 
same binding energy for H, D, H2, HD, and D2. We also verify a 
posteriori that a small difference of 50 K between these binding 
energies has a very limited influence on the deuteration of water. 

To take the binding energy distribution of light adsorbed 
species on ices into account, we assume the binding energy of 
H, D, H2, HD, and D2 relative to ASW E^asw as a f ree param- 
eter between 400 and 600 K. We consider constant binding en- 
ergies of heavier species, following several experiments. Table|4] 
lists the binding energies of selected species. We assume that the 



deuterated species have the same binding energy as their main 
isotopologue. 

Table 4: List of selected species and binding energies relative to 
amorphous solid water ice. 



Species 


E h (K) 


H 


400 - 600 " 


H 2 


400 - 600 a 


C 


800* 


N 


800* 


O 


800* 


CO 


1150 c 


co 2 


2690 d 


2 


1000 c 


3 


1800 e 


N 2 


1000 c 


CH 4 


1300 f 


NH 3 


1300 f 


OH 


2820 - ; 


H 2 


5640 * 


H 2 2 


5640 s 


H,CO 


2050 f 


CH 3 OH 


5530 c 



Refer ence s. {a) Hornekaer etal.1 d2005h; [Al-Halabi & van DishoecM 
j2007h: ( *» iTielens & Allamandolal il987h: (c) ICollings et all d2004l) : 
(rf) ISandford & Allamandolal dl990l>; M : ICuppen & Herbstl d2007h ; (/) 
iGarrod & Herbstl d2006h : M ISpeedv et alJ dl996h 

H2 is, by about four orders of magnitude, the most abundant 
molecule in molecular clouds. Most particles that accrete onto 
the surface are, therefore, H2 molecules. At low temperatures, 
H2 would become the most abundant icy molecule if binding 
energies were computed relative to a water ice subs t rate a lone. 
However, microscop ic models by [Cuppen & Herbstl d2007h and 
ICuppen et all {2009) have shown that the total binding energy of 
a adsorbate relative to a substrate is given by the additive energy 
contribution of the occupied neighbouring sites. Therefore, an 
H2 Q:H2 mixture must be considered for computing an effective 
£ h . IVidari etaD (Il99lh gave an estimate of the binding energy 
of H on an H2 i ce, w hich is about 45 K at 10 K. Following 
IGarrod & Paulvl (1201 ll) . we compute the effective binding en- 
ergy of each species ; from the fractional coverage of H2 on the 
surface, P(H2) 

E b (i) = (1 - P(U 2 )) ■ £„, wat (0 + P(H 2 ) ■ E b H2 (i) (14) 

where E bwa i(i) is the binding energy of the species i relative 
to water ice, the values of selected species are listed in Table 
[4] To deduce E h j_[ 2 of species other than H, we apply the same 
scaling factor as for atomic hydrogen (E b Y[ 2 (H)/E b ASW(H))- 
The increase in the H2 abundance in the mantle tends to decrease 
the effective binding energy Eb(i) of physisorbed species i. 

2. 7. Physical model 

As described in the introduction, w ater ice is formed a t visua l 
extinctions A v higher than 3 mag dWhittet et al.l 1 19881 1200 lb . 
whilst gas pha se CO is detected ab ove a visual extinction thresh- 
old of 2 mag (Frerking et al. 1982) and PDR models show that 
hydrogen is already mainly molecular at Ay < 1 mag. Water 
ice is, therefore, thought to be formed when the gas is already 
molecular. 

We consider a two-step model. For each set of input param- 
eters, we first compute the abundances of gas phase species as- 
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suming the element abundances shown in Tableland consider- 
ing the gas phase network described in section l2~2l For this, we 
assume the steady state values, which are reached after 5 x 10 6 
yr when the density is 10 3 cirT 3 and after 10 6 yr at 10 4 cm -3 . 
This first step is meant to describe the molecular cloud and the 
initial pre-collapse phase before the formation of the ice bulk. 
In practice, we assume that the timescale to reach the chemical 
equilibrium is shorter than the dynamical timescale for the gas to 
reach the prestellar core conditions. These gas phase abundances 
are then considered as the initial abundances for the gas-grain 
modelling. In the second step, we allow gas and grain surface 
chemistry to evolve whilst physical conditions (density, temper- 
ature, visual extinction) remain constant. 



2.8. Multi-parameter approach 

Following our previous works (TCK12a, TCK12b), we consider 
a multi-parameter approach by considering free physical, chem- 
ical, and surface input parameters. This approach allows us to 
study their influence on the formation and the deuteration of key 
species. 

The total density of H nuclei in molecular clouds, where 
interstellar ices are thought to be formed, show typical varia- 
tions from 10 3 cm 3 at the edge of clouds to 10 6 cm 3 in denser 
prestellar cores. Temperatures of the gas and grains also show 
variations depending on the location inside the cloud. Here, we 
consider three fixed temperatures, by assuming that the gas and 
grain surface temperatures are equal, T g - Tj. We study the in- 
fluence of the visual extinction on the formation and the deutera- 
tion of ices by considering A y as a free parameter. As previously 
explained in section l2~2l the ortho/para ratio of H2 is considered 
as constant. 



Table 5: List of the input parameters and the values range ex- 
plored in this work. Bold values mark the values adopted in the 
reference models (see text). 



In section 12.61 and in TCK12a and TCK12b, we showed 
the importance of considering distributions of several grain sur- 
face parameters (grain diameter ad, binding energy Et, of light 
species, diffusion to binding energy ratio E c //Eb). We therefore 
consider them as free parameters, the range of values are listed 
in Table |5] 

We keep fixed the following other parameters: 

- distance betw een two sites d x = 3.1 A , corresponding to a high- 
density ASW dJenniskens et al.lll995h ; 

- cosmic -ray ionization rate ( = 3x 10~ 17 s _1 ; 

- interstellar radiation field (ISRF) Fjsrf = lxlO 8 photons cm -2 
s" 1 ; 

- cosmic ray induced radiation field (CRH2RF) Fcrhirf = 
1 x 10 4 photons cirT 2 s _1 ; 

- scaling factor in multiples of the local interstellar field Go = 20, 
which represents an average value between low interstellar radi- 
ation fields seen in m olecular clouds, such as B68 (Go = 0.25-1, 
iBergin & Sne1iF2 002). and high ISRFs see,n for e xample, in star- 
formi ng molecular clouds in Orion (Go = 10 4 , iGiannini et al.l 
l2000h . 

In the model grid, some of the parameter values are inconsis- 
tent with each other and do not necessarily reflect realistic phys- 
ical models. We use this model grid in order to systematically 
study the influence of each parameter on the formation and the 
deuteration of interstellar ices. For this purpose, we compute the 
mean values and 1-sigma standard deviations of species abun- 
dances and deuterations either in the gas phase or grain surfaces, 
following the method by Wa kelam et al.l (120 1 Oh . Table[5]lists all 
the free parameters and their ranges of values explored in this 
work. 



Input parameters 



Values 



Physical conditions 

T g = T d 
A v 

<r 

ISRF 

Go 

Grain surface parameters 

a d 

E b m 

Ed I Ei, 
d, 

Chemical parameters 

H 2 o/p ratio 



10 3 - 10 4 - 10 5 - 10 6 cm" 3 
10 - 15 - 20 K 
- 1 - 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 mag 

3 x 10~ 17 s _I 
1 x 10 8 photons cm -2 s" 1 
20 

0.1 -0.2 -0.3 /mi 
400 - 500 - 600 K 
0.5 - 0.65 - 0.8 

3.1 A 

3 x 10- 6 - 3 x 10- 4 - 3 x 10- 2 - 3 



3. Results 

In this section, we present the results in two steps. First, we con- 
sider the formation of water and other major ice species, and 
second, we focus on the deuteration. 

Specifically, in section |3~T1 we validate our model. For this 
purpose, i) we compare our predicted gas phase abundances as 
function of the visual extinction with PDR model predictions; ii) 
we compare the predicted abundances of gaseous and solid water 
as function of the visual extinction with published observations; 
iii) we discuss the multilayer formation of interstellar ices for 
a set of three reference physical conditions. We emphasize that 
our approach does not pretend to describe the whole evolution of 
the cloud. The reference models are meant to quantify the influ- 
ence of specific physical conditions, which are likely to describe 
different evolutionary stages, on the ice chemistry. 

Second, after validating the model, we focus on the water 
deuteration. In section I3~2l we emphasize the importance of the 
CO depletion on the deuterium fractionation of the reference 
models. Then, in section l3~3l we perform a multiparameter study 
that allows us to evaluate the influence of each model parameter 
on the deuterium fractionation of water. Comparing the model 
predictions with the observations allows us to constrain a range 
of values for the chemical and physical parameters. 

3. 1. Validation of the model 



3.1 .1 . Initial gas phase abundances 



As described in section I2771 we adopt a two-step model. The ice 
formation phase is followed by considering initial abundances 
that are computed from steady-state calculations of gas phase 
chemistry. These initial abundances depend on the density and 
the temperature, but mainly on the visual extinction Ay (here- 
after, Ay means the edge-to-centre visual extinction, half of the 
observed visual extinction), because photolytic processes play a 
significant role at low visual extinctions. In fact, the abundance 
of most gas phase species weakly depends on the density and the 
temperatures (in the range of values considered in this work) but 
are essentially a function of the visual extinction. Therefore, we 
compare the gas-phase abundances of key molecules for ice for- 
mation (H, D, C, O, CO, O?) betw een our model and the Meudon 
PDR code dLe Petit et ail l2006h for a total density n H = 10 4 
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Figuref3]shows the depth-dependent chemical abundances of 
gas phase H, D, O, C, O2, and CO, for riff = 10 4 cirT 3 computed 
by the PDR code in the low visual extinction regime. These 
abundance profiles are in relative good agreement with the initial 
abundances computed with our model. The abundances differ by 
20% for most abundant species and within one order of mag- 
nitude when species display low abundances (C and O2). This 
validates the use of the H2, HD, and CO self-shieldings in our 
code. The other differences are due to the different chemical net- 
works and physical conditions between the two codes (the PDR 
and ours). The decrease in UV flux with increasing Ay increases 
the abundance of molecules such as H2, HD, and CO, and de- 
creases the abundances of H, D, and then C, and O. D abundance 
is governed by the high photodissociation of HD at low Ay and 
by the formation of H\ isotopologues at higher Ay, and shows 
an abundance minimum at Ay = 1 mag. 

Molecular oxygen is formed at higher visual extinction, and 
it reaches its maximal abundance at Ay > 5 mag. Our chemical 
network predicts a high O2/O abundance ratio (~ 1 - 1 .5). In fact, 
the O2 abundance is highly uncertain because it depends on 1) 
the rate coefficients of some key reactions displaying high uncer- 
tainti es, and 2) the elemental abundances of car bon and oxygen 
(see IWakelam et al.ll20Tot iHincelin et alJl201 lh . Moreover, ob- 
servations carried out with the SWASS, Odin, and Herschel space 
telescopes have revealed th at O2 is not abundant in molecular 
cloud s (X(0 2 ) < 10" 7 , see Gold smith et"alll201 it ILiseau et all 
|2012|) . Modelling the formation of ices with high O2 abundances 
is not necessarily realistic, and we, therefore, investigate the in- 
fluence of the gas phase O2/O abundance ratio on the formation 
of ices in section U. 1.31 




1 2 3 4 5 

Visual extinction (mag) 



Fig. 3: Initial abundances of C, O, CO, and O2 as function of 
the visual extinction computed with the Meudon PDR code at 
n H = 10 4 cm" 3 , Go = 20. 



Concluding remarks. As shown in Fig. [3] the chemistry in the 
gas phase is known to strongly depend on the visual extinction 
Ay. The steady-state gas phase abundances computed with our 
model are in good agreement with the results of PDR simula- 
tions. The use of a two-step model might overestimate the abun- 
dance of gaseous O2 affecting ice formation. Therefore, the in- 
fluence of the O2/O ratio on ice formation needs to be investi- 
gated. 



3.1 .2. Depth-dependent water abundances 

Figure [4] shows the final (at 10 7 yr) abundances of water ice and 
vapour as functions of visual extinction Ay for physical condi- 
tion representative of molecular clouds (e.g. = 10 4 cm -3 , 
T g = Td — 15 K). The final abundances can be divided into two 
zones: 

i) a photon-dominated layer (Ay < 2 mag) where gas phase 
molecules are photodissociated and ices are photodesorbed. At 
very low visual extinctions, interstellar ices are efficiently pho- 
todesorbed and show a low abundance of ~ 10~ 7 relative to H 
nuclei corresponding to less than one monolayer. The abundance 
of water vapour, given by the balance between its formation in 
gas phase, its photodesorption from interstellar grains, and its 
photo dissociation, reaches s teady state values of 10 s — 10 7 (see 
also lHollenbach et alJl2009h . 

ii) a darker region (Ay > 2 mag) where the decrease in UV 
flux allows the formation of interstellar ices, mainly composed 
of water. At Ay = 2-4, water ice is mainly formed from the 
accretion of O atoms via the barrierless reaction (|2)- Most of the 
oxygen reservoir not trapped in CO is easily converted into wa- 
ter ice, reaching abundances up to 2 x 10~ 4 . The increase in gas- 
phase O2 initial abundance with Ay slightly decreases the final 
abundance of water ice. Indeed, water ice is also formed via the 
accretion of O2 including the formation of hydrogen peroxide. 
These formation pathways involve reactions possessing signifi- 
cant activation energies and low transmission probabilities (see 
Table[2]). Finally, water formation also depends on the grain sur- 
face parameters Ed/Eb and Et,(H). A highly porous case (high 
Ed I Eb ratio combined with a high binding energy of H) strongly 
decreases the diffusion rate of mobile H species, decreasing the 
final abundance of water ice by one order of magnitude (~ 10~ 5 ). 

Gas phase abundance of water is mainly governed by the 
balance between photodesorption and accretion. The decrease 
in UV flux with increasing A v decreases water ice photodesorp- 
tion and therefore the final abundance of water vapour from 10 7 
at Ay = 2 to less than 10~ 9 at Ay = 5. We can note that the 
multilayer nature of grain mantles and the use of a wavelength- 
dependent absorption of UV photon s from ices give similar re- 
sults to theoretical PDR studies (e.g. lHollenbach et alj|2 009). 

In spite of the inevitable approximations of our modelling, 
our predictions are in good agreement with observations of water 
ice and water vapour, as shown in Fig. |4] 

Concluding remarks. For typical conditions representative of 
a molecular cloud, we are able to reproduce the observed Ay 
threshold (Ay ~ 1.5 mag, Ay„t, s ~ 3 mag) and the high abun- 
dance of water ice (X ~ 10~ 4 ). The low abundance of water 
vapour observed in molecular clouds is reproduced at higher vi- 
sual extinctions (4-5 mag). 

3.1.3. Reference models 

Very likely, the different ice components and their relative 
deuteration are the result of a long history where the physical 
conditions evolve. Therefore, a model aiming at reproducing the 
whole set of observations should take this evolution into account. 
However, before embarking on such a complicated modelling, it 
is worth while and even important to make clear what the char- 
acteristics (specifically, solid species and deuteration) are at each 
evolutionary step. To this end, we consider three "reference" 
models in the following: i) translucent cloud model, forming 
H2O and CO2 ices; ii) dark cloud model, allowing the forma- 
tion of CO ice; iii) dark core model, showing high depletion of 
CO. The chemical composition of grain mantles for these refer- 
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Fig. 5: Fractional ice abundances of main species formed from O, O2, and CO for three reference models with "normal" gas phase 
abundances computed with the chemical network presented in section I2.2I (left), and for a gas phase [O2MO] abundance ratio 
artificially decreased by a factor of 10 (right): top) = 10 4 cm -3 , T — 15 K, Ay = 2 mag, middle) tin = 1 X 10 4 cirT 3 , T = 10 K, 
Ay = 4 mag, bottom) n# = 10 5 cirT 3 , T = 10 K, Ay — 10 mag. Values of other parameters are bold values of Table|5] 



ence models are shown in Figure |5]for the two sets of runs with 
different [Chl/tO] abundance ratios (see section l2~!2l . 

i) Translucent cloud model where ices are mainly composed 
of H 2 and C0 2 : n H = 10 4 cm 3 , T = 15 K, Ay = 2 mag (cor- 
responding to an observed visual extinction of 4 mag). For the 
conditions considered in this model, grain mantles are mainly 
composed of water and carbon dioxide. After H 2 , the most abun- 
dant species that accrete onto grain surfaces are H, O, and CO. 



For these conditions, the initial O2/O abundance ratio is low 
(5 x 10~ 4 ). The decrease in the O2/O ratio, therefore, does not 
affect the ice composition significantly. 2 is involved in the for- 
mation of O and CO in the gas phase, and its artificial decrease 
favours O instead of CO and increases the formation of water ice 
by 25%, increasing the ice thickness. The relatively high grain 
temperature (15 K) allows accreted particles to diffuse efficiently 
and form H2O and CO2. At this visual extinction, the UV flux 
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Fig. 4: Final abundances (10 7 yr) of water ice (H20 /ce , red) 
and water vapour (H20 gas , blue) as function of the visual ex- 
tinction A Y for typical molecular cloud conditions («# = 10 4 
cm -3 , T = 15 K) and including the variation of all other in- 
put parameters. H a tched boxes refer to water ice observed by 
iPontoppidan et alj d2004l) tow ards the Serp e ns SM M4 regions 
and water vapour observed by ICaselli et al.l (1201 Ob towards the 
L1544 prestellar core (the upper limit represents water abun- 
dance in the external part whilst the lower value represents water 
abundance in the central region). 



irradiating grain surfaces is high, and volatile species, such as 
atoms or light stable molecules (CO), are efficiently photodes- 
orbed. 

Water is the main ice component because it is formed from 
the banierless reactions ([T} and (fj} whilst its other reaction 
routes are negligible (less than 0. 1 %). Carbon dioxide is mainly 
formed from the hydrogenation of the O...CO complex and not 
by the direct reaction between CO and OH. Indeed, the grain 
temperature is not high enough to allow a high diffusion of CO 
and OH, because of their high binding energy. Instead, O atoms 
that are less attractively bound can meet CO atoms to form the 
O...CO van der Waals complex. In turn, O...CO readily reacts 
with H to form a hot HO. .CO* complex that forms C0 2 + H. 
The accretion of CO and O also allow the formation of formalde- 
hyde and methanol, but only in low abundances, lower than 1 % 
compared to solid water whilst a low fraction of hydrogen per- 
oxide is predicted for a "normal" O2/O ratio. Indeed, they are 
formed by reactions involving either heavy atoms or significant 
activation barriers. 

The overall abundance of CO2 relative to H2O decreases 
with time, varying from 30% at the beginning to less than 1% 
at the end. The decrease in the CO2 abundance is due to the 
increase in H abundance once interstellar ices start to form, be- 
cause of the high photodesorption. 

Although H2O and CO2 are formed in high abundances, O 
and CO are not totally depleted on grains. They still show high 
gas-phase abundances (about 5 x 10 5 relative to H nuclei) at the 
end of the ice formation. The predicted CO2 abundance relative 
to water ice is lower than the observations. However, CO2 for- 
mation is very sensitive to the grain temperature (governing the 
diffusion of O atoms) and to the hydrogen abundance (governing 
the reaction rate between O and H atoms). A higher temperature, 
higher visual extinction, and/or higher density would, therefore, 
tend to increase the abundance of CO2 relative to water. 



ii) Dark cloud model where CO ice starts to form with water, 
forming a CO:H 2 mixture: n H = 1 x 10 4 cm 4 , T = 10 K, 
Ay — 4 mag (corresponding to an observed visual extinction of 
8 mag), t = 10 6 yr. The low temperature limits the diffusion of 

atoms and heavier species whilst the higher visual extinction 
allows CO molecules to stay bound on grain surfaces even if 
they do not react to form CO 2. Therefore, the main components 
of grain mantles are H2O, mainly formed via the reaction (f2}, 
and CO. 

CO shows an overall abundance relative to water of about 
80 %, which is more than two t imes higher than observed abun- 
dance ratios of ices (10-30 %, IWhittet et ai]|2007t IQberg et al.l 
2011). This is due to the high O2/O abundance ratio (~ 0.5 
throughout the calculation) that allows the formation of H2O2, 
via banierless hydrogenation reactions, instead of water. H2O2 
is easily trapped within grain mantles before reacting because 
the reaction destroying hydrogen peroxide has a low transmis- 
sion probability (see Table [2). The decrease in the O2/O ratio by 
one order of magnitude increases the formation of water, giv- 
ing a CO/H2O abundance ratio of 30%, whilst the H2O2 abun- 
dance decreases by one order of magnitude. Hy d rogen peroxide 
has recently been observed by Bergman e t al.l (1201 lb towards 
the p Oph A dark cloud, confirming its formation in ices and 
its subsequent sublimation via non-photolytic processes. The 
low grain temperature favours the formation of formaldehyde 
and methanol via CO hydrogenation, compared to the formation 
of CO2. Formaldehyde shows an unexpected abundance higher 
than 8 % compared to water ice. Again, the formation efficiency 
of formaldehyde and methanol strongly depends on the grain 
temperature and on the total density. For example, a higher tem- 
perature and/or higher density would favour CO2 formation in- 
stead of H 2 CO and CH3OH. 

A visual extinction of 4 mag decreases the final gas-phase 
abundances of O and CO, compared to Ay = 2 mag. Indeed, gas- 
phase CO reaches an abundance of 5 x 10~ 6 whilst O abundance 
decreases to 10~ 7 at 10 7 yr. The difference in the two abundances 
is due to the efficient destruction of O atoms forming solid wa- 
ter, whilst most CO molecules do not react and are still able to 
photodesorb. 

iii) Dark core model where most of CO is depleted, allow- 
ing the formation of pure CO ice and solid formaldehyde and 
methanol: Hh = 10 5 crrT 3 , T = 10 K, Ay = 10 mag (corre- 
sponding to an observed Ay of 20 mag). The low temperature, 
the high visual extinction, and the high density allow a signifi- 
cant trapping of CO molecules in the inner part of grain mantles 
whilst the formation efficiency of CO2 is low owing to the low 
temperature. H2CO and CH3OH are formed via hydrogenation 
reactions that have high activation barriers. They are, therefore, 
mainly formed in the outer part of the ice when the CO depletion 
is high, which allows an efficient hydrogenation (see TCK12a). 

At this visual extinction, solid water is less abundant than 
hydrogen peroxide because the initial gas-phase O2 abundance 
is higher than the abundance of atomic oxygen (O2/O between 

1 ~ 3 throughout the calculation). Water formation is mainly 
formed through the formation of H2O2, which involves reac- 
tions having high activation barriers. Hydrogen peroxide is most 
likely trapped in the bulk before forming water owing to the rel- 
atively high density. In this case, solid CO reaches an absolute 
abundance relative to H nuclei of 10~ 4 , water and CO2 show 
lower abundances (10~ 5 and 4 x 10~ 6 ), whilst formaldehyde and 
methanol show an abundance of 2 x 10~ 6 and 8 x 10~ 7 , respec- 
tively. Decreasing the O2/O abundance ratio to 0.3 strongly in- 
creases the formation efficiency of water and CO2 instead of 
hydrogen peroxide. For this case, water is almost as abundant 
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as CO. Gas-phase abundances of O, O2, and CO decrease with 
time and show high depletions (with final abundances lower than 
10~ 10 ), whereas H abundance remains constant. 

Concluding remarks. The chemical composition of ices is very 
sensitive to the physical conditions and to the initial abundances. 
Most of the observed ice features are reproduced by the refer- 
ence models: the water-rich ice seen at low visual extinctions is 
also composed of abundant CO2, whilst the abundances of solid 
CO (and H2CO and CH3OH) gradually increase with the visual 
extinction and the density. The increase in the O2/O ratio de- 
creases the abundance of water because its formation from O2 
involves reactions that have significant activation barriers (see 
Table 0. However, since it is likely that O2 abundance remains 
low in dark clouds, water formation seems to be efficient in a 
wide range of physical conditions. Therefore, the study of the 
water deuteration needs to include the variation of several phys- 
ical parameters. 



3.2. CO depletion and molecular deuteration 




x 

X 



0.001 - 



1 10 100 

CO depletion 

Fig. 6: Solid HDO/H2O (blue) and gaseous D/H ratio (orange) as 
function of the CO depletion factor f D (CO) - n g (CO)ln gMi (CO), 
which increases with time, for the three reference models: 
translucent cloud region (dotted), dark cloud region (dashed), 
dark core region (solid). 



As discussed in TCK12b, the deuteration of solid species 
strongly depends on the values of the CO depletion and the den- 
sity at the moment of their formation. The initial densities of H 
and D are roughly constant regardless of the total density. Their 
abundance relative to H nuclei decreases with increasing At 
low densities, the increase in D abundance is limited by the weak 
deuterium reservoir, whilst the D/H ratio is able to strongly in- 
crease at higher densities. The gas and grain temperatures, as 
well as the visual extinction, also affect the evolution of the D/H 
ratio because they influence the desorption rate of H and D. In 
summary, since all these quantities vary with time, deuteration 
is not necessarily constant within grain mantles. 

The influence of the CO depletion, the density, the tempera- 
tures, and the visual extinction on the water deuteration is shown 
in Figure[6]for the three "reference" models described in the pre- 
vious section. At the beginning of the gas-grain calculation, CO 
has already reached its maximal abundance (~ 1CT 4 , see Fig. [3]). 
With time, the CO molecules freeze-out onto grains, decreasing 
the gas-phase CO abundance and increasing the CO depletion 
factor. The decrease in the temperatures from 15 to 10 K and the 
increase in the visual extinction from 2 to 4 mag increase the gas- 
phase D/H ratio from 0. 1 % to 0.4 % at f D (CO) = 1 . The increase 
in total density from 10 4 to 10 5 crrT 3 increases the final D/H ra- 
tio from 1 % to 10 %. The decrease in the gas-phase O2/O abun- 
dance ratio does not modify the water deuteration for the two 
low-density cases. However, it slightly decreases the HDO/H2O 
ratio for the dense core model by a factor of two because water 
is more efficiently formed at low CO depletion (see Fig. O when 
the D/H ratio is low. 

Concluding remarks. As previously noted, water deuteration is 
largely influenced by the gas D/H atomic ratio. Consequently, 
the largest deuteration is obtained where the CO depletion and 
the gas-phase D/H ratio are high, namely in the latest and less 
efficient phases of water formation, represented by the reference 
model iii). Comparison with the observed deuterium fractiona- 
tion towards IRAS 16293 shows that a part of water ice should 
have formed during a dark and/or dense phase. In the following 
section, we investigate the impact of key parameters on water 
deuteration. 



3.3. Physical/chemical parameters and water deuteration 

In this section, we study the influence on ice deuteration of sev- 
eral physical and chemical parameters which play a key role. 
Each figure presented in this section shows the influence of one 
(or two) parameter(s) at a time. For each value of the studied pa- 
rameters, mean value and standard deviation of absolute abun- 
dances and deuterations induced by the variation in other input 
parameters are computed. Comparisons between the evolution 
of the mean deuteration induced by the variation in the studied 
parameter and the standard deviation caused by other parameters 
allow us to deduce the importance of that parameter on water ice 
deuteration. 

The predicted deuteration levels of water ice are compared 
with the observed HDO/H 9O and D7O/H7O ra tios towards the 
hot corino of IRAS 16293 Icoutens et al .1120121) to constrain the 
input parameters that reproduce the observations best. For this 
purpose, we assume that the entire bulk of interstellar ices des- 
orbed in the hot corino and the observed deuteration reflects the 
deuteration in ices (see Introduction). 

3.3.1 . Influence of the H 2 o/p ratio 

As described in section 2.2, the H2 opr influen ces the deutera- 
tion o f gas-phase sp ecies (including ato mic D) JWalmsle v et al.l 
2004; iFlower et al.ll2006t iPagani et al.ll2009l) . thereby affecting 
the deuteration of solid water. Figure |7] shows the HDO/H2O 
and D2O/H2O ratios in grain mantles as a function of time for 
four values of H2 opr, including the variation in all other param- 
eters, except the visual extinction range, which is limited to 2 - 
10 mag (where most of water is believed to form, see Fig. |4). As 
anticipated, the H2 opr governs the deuteration of water ice via 
its influence on the abundance of H3 isotopologues and atomic 
D, for values higher than about 3 x 10~ 4 . An increase in H2 opr 
between 3 x 10" 4 and 3 decreases the HDO/H 2 and D 2 0/H 2 
ratios by 2.5 orders and 5.5 orders of magnitude, respectively. 
Furthermore, the ortho/para ratio of H2 is a key, even the most 
important, parameter for water deuteration. The decrease in the 
deuteration due to the increase in the H2 opr is much greater than 
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the standard deviation induced by the variation in other parame- 
ters. 




10 s 10 6 
Time (years) 

Fig. 7: Deuteration of water ice for four values of H2 ortho/para 
ratios: 3 x 1(T 6 (red), 3 x 1(T 4 (blue), 3 x 1(T 2 (green), and 
3 (yellow) including the variation in all other input parame- 
ters, except the visual extinction range that is limited to 2 - 
10 mag. Hatched bo xes refer to water deuteration observed by 
ICoutens et all (120121) towards IRAS 1 6293 . 



Comparisons with the observations by ICoutens et al.l (1201 2|) 
clearly suggest that a low opr of H2, lower than 3 x 10~ 4 , is 
needed to reproduce the observed deuterium fractionation. 

3.3.2. Influence of the total density 

As discussed in section 13.21 the gas-phase D/H ratio is a func- 
tion of the total density njj because high densities allow the 
gas-phase D/H ratio to increase with the CO depletion. Figure 
[8] shows the deuteration of HDO and D2O for the four consid- 
ered densities, using an H2 opr of 3 x 10~ 6 , two temperatures (10 
and 20 K), and a high visual extinction (10 mag), including the 
variation in grain surface parameters. 

At low densities (n# < 10 4 cm -3 ), water deuteration is con- 
stant and low with time because a significant part of HD is al- 
ready trapped in atomic D, before the CO depletion. At higher 
densities {tin > 10 5 cm 4 ), the efficient increase in the gaseous 
D/H ratio allows an increase in the deuteration of water up to 
10%. In these cases, HDO is mostly located in the outer part 



of grain mantles, and D2O is only located in the outermost lay- 
ers, whilst H2O is present throughout the mantle bulk. For this 
reason, the HDO/H2O is limited and cannot reach the final gas- 
phase D/H ratio (up to 50 % at n H = 10 6 cm 4 ). 

In summary, the total density plays a key role in the deuter- 
ation of water. Indeed, the variation in the total density njj be- 
tween 10 3 and 10 6 cm 3 influences the HDO/H 2 and D 2 /H 2 
ratios by 2 and 3.5 orders of magnitude, respectively. 

The two observed deuteration ratios can be predicted with to- 
tal densities njj between 10 4 and 10 5 cm 4 at 10 K and between 
10 s and 10 6 cm 4 at 20 K, regardless of other grain surface pa- 
rameters. In the following section, we study the case h# = 10 4 
cm -3 in more detail to investigate the effect of the temperature 
and the visual extinction on water deuteration. 



3.3.3. Influence of the temperature and visual extinction 

Temperatures either in the gas-phase and on grain surfaces also 
influence the deuteration of water. First, an increase in the gas- 
phase temperature enhances the reactivity of endothermic reac- 
tions, which can hydrogenate back Hi isotopologues, decreas- 
ing the abundance of gaseous D atoms. Second, an increase 
in the grain temperature severely increases the desorption rate 
of volatile species, such as atomic H or D. The abundance of 
gaseous D is limited by the low deuterium reservoir, whereas 
the abundance of atomic hydrogen can increase up to two orders 
of magnitude, decreasing the gaseous D/H ratio. Water is mainly 
produced via reactions involving atomic H and D, increasing gas 
and dust temperatures, hence, decreasing the water deuteration 
on grain surfaces. 

Figure [9] shows this effect by presenting the final HDO/H2O 
and D2O/H2O ratios (at 10 7 yr) with the visual extinction, for 
the three temperatures, considering an H2 opr of 3 x 10~ 6 , a total 
density of 10 4 cm 4 , and including the variation in grain surface 
parameters. The increase in the (gas and grain) temperatures be- 
tween 10 and 20 K slightly decreases the water deuteration ap- 
proximately by a factor of 3 for HDO and of 10 for D2O at visual 
extinctions higher than 2 mag. The evolution of deuteration is the 
same order of magnitude as the standard deviations induced by 
the variation of other parameters. 

The H and CO abundances in the gas-phase increase with de- 
creasing Ay because of the increase in the photodesorption rate, 
decreasing the gaseous D/H ratio. Visual extinction, therefore, 
influences water deuteration at low visual extinctions where wa- 
ter starts to form (2mag < Ay < 5 mag). As can be seen in Fig. [9] 
HDO/H2O and D2O/H2O ratios increase by one and two orders 
of magnitude respectively, between Ay — 2 and 5 mag, . 

Comparison with observations shows that observed HDO 
and D2O fractionations are reproduced for a low temperature of 
10 K and for visual extinctions Ay higher than 4 mag (Av, bs = 8 
mag), implying that deuterated water needs to be formed in dark 
regions if = 10 4 cm" 3 . 



3.3.4. Influence of grain surface parameters 

The grain surface parameters {ad, Ed I Eh, Et,(H)) can also influ- 
ence the formation of interstellar ices. As discussed in TCK12a, 
the absolute abundance of the main ice constituents formed on 
grain surfaces (water, formaldehyde, methanol for instance) de- 
crease with Ed I Eb whilst the grain size does not affect the overall 
abundance (but only the ice thickness). Moreover, it is also seen 
that the absolute abundance of water slightly decreases with the 
binding energy of volatile species. 
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Fig. 8: Deuteration of water ice for four values of total density 10 3 (red), 10 4 (blue), 10 5 (green), 10 6 (yellow) cm -3 , and two 
temperatures: 10 K (left panels), 20 K (right panels), using an H2 opr of 3 x 10~ 6 , a visual extinction of 1 mag , including the 
variation of all other input parameters. Hatched boxes refer to water deuteration observed by ICoutens et al.l (12012) towards IRAS 
16293. 



The water deuteration is slightly affected by the variation of 
the grain surface parameters. As shown in Fig. |9] the standard 
deviations given by the variation of the grain surface parameters 
remain lower than the evolution of the water deuteration lead by 
the variation in the temperature and the visual extinction. The 
water deuteration slightly decreases with the grain size due to 
the decrease in the CO depletion efficiency (since the accretion 
rate is inversely proportional to the grain size), limiting the in- 
crease in the deuterium fractionation. On the other hand, the wa- 
ter deuteration slightly increases with Ej/Eb and Et,(H) because 
the formation efficiency of H2O decreases with the diffusion en- 
ergy of volatile species, increasing the overall fractionation. 

3.3.5. Concluding remarks 

Comparing the average value of the deuterium fractionation with 
the standard deviation induced by the variation of other param- 
eters allows us to claim that the H2 opr is the most important 
parameter for water deuteration, followed by the total density. 
The visual extinction and the temperature also influence water 
deuteration but more weakly. 

The comparison of our model predictions with observations 
allows us to constrain some parameters involved in the formation 
of deuterated water seen towards IRAS 16293. The observations 
are reproduced for 



- a H2 opr lower than 3 x 10~ 4 ; 

- a density lower than 10 5 crrT 3 if T = 10 K and Ay > 4 mag; 

- a higher density nj_[ between ~ 5 x 10 4 cirT 3 and 10 6 crrT 3 if 
T = 20 K. 

4. Comparisons with previous models 

In this section, we compare our model predictions for the gas- 
phase deuteration and for water ice deuteration with previous 
models. The atomic D/H ratio in the gas-phase depends on the 
deuteration of H^ . W e compare our pred ictions from steady state 
models with those bv lFlower et al.l (120061) . specifically their Fig. 

5. We obtain the same influence of the H2 opr on the deuteration 
of Hi , with similar fractionation values. 

To our knowledge, the only comprehensive theoretical study 
focused on the deuterated water fo rmation on grain surfaces was 
performed by Caza ux et al 1 (1201 11) . even if other works also in- 
clude d the deuterati o n of water ice in thei r astrochemical model 
(i.e. iTieleni 119831: ICaselli et al.l 120021: IStantcheva & Herbstl 
120031: iBell et alJl201 lh . lCazaux et alJ d201 ll) considered a static 
stage followed by a free-fall collapse phase. However, in their 
model, the deuterated water is practically formed only during 
the first static phase. Thus, it is worth comparing their predic- 
tions with ours. 

The main difference is that they found that the water deuter- 
ation is highly temperature-dependent and lower than our pre- 
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Fig. 9: Final deuteration of water ice (at t > 10 yr: representa- 
tive of a typical molecular cloud age) as function of visual ex- 
tinction for three values of temperatures: 10 (red), 15 (blue), 20 
(green) K, using a total density of 10 4 cirT 3 , an H2 opr of 3x 10~ 6 
and including the variation of all other input parameters. The 
vertical line refers to the Ay threshold of water ice observed by 
Whit tet et al. ( 19881). Hatched bo xes refer to water deuteration 
observed bv lCoutens etall d2012l) towards IRAS 16293. 



dictions at low ( < 15 K) temperatu res (Fig. [8] of this work 
versus Fig. 3.c of ICazaux et al. This is due to a differ- 

ent approach in the reaction probability computation and the 
use of the endotherm ic O+H2 reaction in the chemical network. 
ICazaux et alj {2011) assumed a competition between the reac- 
tion and the diffusion. In their model, if the diffusion timescale 
of the reactants is larger than the transmission timescale of the 
reaction, the reaction occurs regardless of its activation barrier. 
Thus, at low temperatures, the mobility is low and the reaction 
always occurs. In contrast, in our work, the reaction rate is given 
by the product of the collision rate of the two reactants and the 
reaction probability P r (see Sec. 2.1 and 2.4), giving a much 
smoother dep endence on the temperature. As a consequence, in 
ICazaux et al.l d201 lb . H2O is formed, at low temperatures, via re- 
actions involving H2 because of its large abundance, whilst HDO 
is formed via reactions involving atomic D. The solid HDO/H2O 
ratio, therefore, scales with the gaseous [D]/[H2] ratio (~ 10~ 5 ). 
In our work, water is always mostly formed via the O+H and 
OH+H reactions, and the deuteration scales with the [D]/[H] ra- 
tio (up to 10 _1 ). At higher temperatures (> 15 K), both mod- 
els agree since water is formed via the same reactions, namely 



O+H and OH+H (and their deuterated counterparts). This com- 
parison shows that the chemical route for the water formation is 
the major actor in the game, whilst other differences in the two 
models, including the physical evolution treatment, play a minor 
role in the water deuteration. However, we should point out that 
the O+H2 reaction cannot occur at low temperatures given its 
endothermicity. Therefore, it is unlikely that water deuteration 
shows a such strong temperature dependence. 

5. Discussion 

In the previous paragraphs, we have shown that the deuteration 
of water and other molecules strongly depends on the H2 opr and 
the total density, but also on the visual extinction and the temper- 
ature where solid molecules are formed. The high deuteration of 
water (HDO/H 2 > 1 % and D 2 0/H 2 > 0.01 %) observed to- 
wards the low-mass protostar IRAS 16293 (ICoutens et al.ll2012h 
can only be reproduced with an H2 opr that is lower than 3 x 10" 4 , 
suggesting that most of deuterated water is formed at low H2 
opr. The ortho-para ratio reaches its steady-state value on a 
timescale proportio nal to the total de nsity, of about 10 7 yr at 
n H = 10 4 cm" 3 dFloweretal .1120061) . If the H 2 molecule has 
an initial opr value of 3 upon its formation on grain surfaces, 
the high deuteration observed around IRAS 16293 would sug- 
gest that water observed in this envelope has been formed in an 
"old" molecular gas, i.e. a gas old enough to show a low H2 
opr at the moment of the formation of interstellar ices. However, 
it is possible t hat the ortho-to-p ara conversion also occurs on 
grain surfaces dLe B ourlot 2000) b ut with an uncer t ain rate (see 
Sugi moto & Fukutanill201 U IChehrouri et al.ll201 U lHama et al.l 
2012). The decrease in the H2 opr could occur, therefore, faster 
an d the molecular clou d age of 10 7 yr deduced from the estimate 
bv lFlower et al.l (120061) should only be used as an upper limit. 

The observed deuterium fractionation of water is only repro- 
duced for dark conditions (high visual extinctions and low tem- 
peratures if iih = 10 4 cirT 3 or higher densities). Therefore, al- 
though IR observations of ices show that water ice starts to form 
at low visual extinctions, the deuterated water observed in IRAS 
16293 should instead be formed in darker regions. Consequently, 
H2O would be present in the inner part of ice mantles, whilst 
most of HDO and D2O molecules should be located in the outer 
layers. A physical evolution, modelling the accumulation of mat- 
ter from diffuse molecular clouds to dense cores, would allow us 
to directly confirm this result. The high observed deuteration of 
water also confirms that water is mostly formed from reactions 
involving atomic H and D and not by molecular hydrogen (see 
Secgji. 

Observations of water, formaldehyde, and methanol vapours 
towards low-mass protostars show that these species have 
different deuterium fractionations. Figure ITOl overplots the ob- 
served deuteration towards IRAS 16293 with mean theoretical 
deuteration values reached at t — 3 x 10 5 yr (i.e., the upper 
limit of the age of prestellar cores, see iBergin&TafaTial 12007) 
as function of the density, for two temperatures, a low H2 opr of 
3 X 10~ 6 and only considering dark regions (A v = 10 mag). The 
comparison of our predictions with observations shows that 

i) our model reproduces the observed HDO/H2O and D2O/H2O 
ratios for densities between 1 x 10 4 and 1 x 10 5 cirT 3 at 10 K 
and between and 2 x 10 4 and 3 x 10 5 cirT 3 at 20 K. Therefore, 
water deuteration is reproduced within a wide range of physical 
conditions representative of molecular clouds but not in the 
too-translucent cloud phase with too low density, too low visual 
extinctions and too high temperatures. 

ii) the observed HDCO/H2CO can be reproduced at higher 
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10 cm ) and lower temperatures (~ 10 



densities (3 x 10 s 
K) as seen in the central regions of prestellar cores. The 
D2CO/H2CO ratio is reproduced at densities higher than 10 6 

cm" 3 (at ~ 5 x 10 6 cm" 3 , TCK12b). 

iii) methanol deuteration (CH 2 DOH/CH 3 OH and 
CHD2OH/CH3OH ratios) proceeds in the outer parts of 
grain mantles when the prestellar core condensation reached 
high densities (> 5 x 10 5 cur 3 ) and low temperatures (10 K). 
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Fig. 10: Averaged deuterations of water (green diamonds), 
formaldehyde (blue triangles), and methanol (red squares) at 
3 x 10 5 yr (representative of a typical prestellar core age) as 
function of the total density (abscissa) and temperature (10 K: 
solid, 20 K: dashed) considering an H2 opr of 3 x 10~ 6 and 
Ay — 10 mag and considering grain s urface paramet e rs as f ree. 
Solid boxes refer to observations by ICoutens et al.l (1201 2b to- 
wards IRAS 16293 for water, and bv iParise et al.l (120061) (and 
references therein) towards a sample of low-mass protostars for 
formaldehyde and methanol. 



In summary, the difference in the deuterium fractionation 
of molecules seen in IRAS 16293 can be explained by differ- 
ent periods of formation. Water is mainly formed first in re- 
gions showing intermediate densities representative of molecular 
clouds (rifj — 10 4 - 10 5 cm -3 where the CO depletion is limited. 
Formaldehyde and methanol are formed subsequently, at higher 
CO depletions when the gas is denser and colder. These predic- 
tions are in good agreement with infrared observations of inter- 
stellar ices presented in the introduction which also show that 
water ice is formed first, at low visual extinctions, whilst solid 



methanol forms late r on at visual extinctions higher than 15 mag 
(IWhittet et al.ll2011l) . 

More generally, the deuterium fractionation of any species, 
believed to be formed mainly in ices, can be used as a tracer 
to estimate the physical conditions at the moment of its forma- 
tion. As shown in Fig. [10] the D/H ratio increases with increas- 
ing density and with decreasing temperature. A low deuteration, 
similar to the deuteration of water, implies that the solid species 
is mainly formed in regions with relatively low densities, even 
if its deuterated isotopologues could be formed mainly in a later 
phase. In contrast, a high deuteration level suggests that the solid 
species is mostly formed in the dense and cold phase where 
the deuteration in the gas-phase (and in particular the atomic 
D/H ratio) is high. Extremely high deuteration, such as methanol 
deuteration, would also imply the existence of abstraction reac- 
tions or other processes enhancing the formation of deuterated 
isotopologues with respect to the main isotopologue (TCK12b). 
H/D exchange reactions between solid methanol and water, for 
exampl e, can also lead to a s elective deuteration of functional 
groups (Ratajczaket al. 2009). 



6. Conclusions 

We have presented a comprehensive study of the formation and 
the deuteration of interstellar water ice carried out with our astro- 
chemical model GRAINOBLE. In addition to the multilayer for- 
mation of grain mantles presented in a first paper (TCK12a) and 
the use of abstraction reactions for the deuteration of formalde- 
hyde and methanol (TCK12b), we introduced a better treatment 
for computing the transmission probabilities of surface reactions 
based on the Eckart model. We also considered wavelength- 
dependent UV photodesorption of ices following molecular dy- 
namics (MD) simulations and new experimental works. 
The main results of this work are the following: 

1) Our model reproduces the abund ance of water ice and the 
visual extinction threshold observed bv lWhittet et al.l (Il988l) . as 
well as the different ice components seen in different molecular 
clouds regions. 

2) Water deuteration strongly depends on the ortho-to-para 
ratio of H2 and the total density, but also, even though more 
weakly, on the gas and grain temperatures and on the visual ex- 
tinction. 

3) The deuteration of water observed towards the low-mass 
protostar IRAS 16293 can only be reproduced by considering a 
H2 opr lower than 3 x 10 and a total density between 8 x 10 3 
and 5 x 10 5 cm 4 . If a low density (10 4 cm" 3 ) is considered, 
a low temperature (10 K) and a visual extinction higher than 
4 mag are necessary. Dark regions (high densities or high visual 
extinctions) are, therefore, needed to reproduce the observations. 

4) Comparison between the observed deuteration of water, 
formaldehyde, and methanol and our theoretical predictions al- 
lowed us to propose the following scenario. Water ice is formed 
first within relatively wide ranges of physical conditions: H2O 
is allowed to start its formation at low densities and low vi- 
sual extinctions, but HDO and D2O are instead formed in darker 
(higher njj and/or Ay) regions. Formaldehyde and then methanol 
are mainly formed subsequently, in the dense and cold prestellar 
cores. 

5) Deuterium fractionation can be used to estimate the values 
of the density and temperature at the moment of formation of 
solid species because of its sensitivity to the physical conditions. 
Low deuteration, similar to water deuteration, implies that the 
species is formed in a wide range of physical conditions, whilst 
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a high deuteration suggests an efficient formation in the centre 
of dense and cold prestellar cores. 

With this work, we have explored the influence of a wide 
range of parameters values on the formation of deuterated ices. 
Now that the influence of the various parameters is clarified, our 
next step will be a more sophisticated model, based on a realistic 
physical evolution of a cloud, in a forecoming paper. 
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Appendix A: Transmission probability 
computations 

A. 1 . Eckart model 

We compute the transmission probabilities of all the re- 
actions involved in the water-producing and t he methanol- 
producin g networks by usi ng the Eckart model (lEckartlll930t 
I Johnston & Heicklenlll962h . In this approach, an approximate 
potential energy surface (PES) is fitted as a function of the zero 
point energies (ZPEs) of the stationary points. The parameters 
needed for computing the transmission probability are the zero- 
point-corrected barrier heights of the forward and reverse reac- 
tions Vf and V r , the frequency of the imaginary mode of the 
transition state vs , and the reduced mass of the reactants ji. 
The Eckart potential can be parametrised as 



U x 



Aexp(^) Bexp(^p) 



1 + 



exp(^) (l + exp(^)) 2 



where 



A - Vf - V r 




1 1 

+ 



/v f vv;j 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



Once one has fitted the potential then the transmission prob- 
ability, P r , may be calculated using 



Pr 



cosh(a- + P) — cosh(a -/J) 
cosh(a + f3) + cosh(<5) 



(A.5) 



where 



An 



1 1 
+ 



8tt 2 



h\v S \ 



VI 



V£ - V f + V r 



6 = 4n 



V f V r 2n 1 



V(^|v S |) 2 
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(A.6) 



(A.7) 



(A.8) 



When the reactants are considered as excited by their formation 
involving an exothermic reaction, E refers to the excess energy 
of this reaction of formation. Otherwise, E is the thermal energy 
of the particles. 

The Eckart model provides a significant improvement over 
square barriers but it can underpredict (or overpredict) the 
transmission probabilities of some reactio ns at low tempera - 
tures compared to more exact methods (see iPeters et al.l 1201 lb . 
However, given the number of surface reactions considered in 
this work, exact quantum chemical computations for all reac- 
tions are not feasible. 



A.2. C0 2 formation 

Based on the experimental work bv lOba et al. who stud- 

ied the formation of CO2 from CO and OH, CO2 is thought to 
be formed via the pathway 



CO + OH trans-HOCO 



trans-HOCO -> cis-HOCO 



cis-HOCO C0 2 + H. 



(A.9) 
(A. 10) 
(A. 11) 



The electronic energy of the intermediate radicals (t-HOCO, c- 
HOCO) and the products (CO2 + H) are lower than for the re- 
actants (CO + OH). Reactions involving t-HOCO and c-HOCO 
posse ss high activation barriers and/or are endothermic (lYu et al.l 
120011) . Therefore, t-HOCO radicals continue to react only if their 
excess energy released by the chemical energy of reaction (1A.9I > 
is sufficient to overcome t he high activatio n barrie rs of reactions 
dATOl l and dATTI As in IGoumans et all (l2008h . HOCO radi- 
cals can also react with H atoms via barrierless reactions to form 
three pairs of products: H 2 + C0 2 , H 2 + CO, or HCOOH. The 
lack of more quantitative data leads us to assume a branching 
ratio of one to three for every product pair. 

In fact, the energy released by reaction (I A. 9) absorbed by 
HOCO radicals can be transferred to the surface before reacting. 
However, the transfer rate of the chemical energy to the surface 
is very uncertain. The absence of formic acid HCOOH and the 
lo w abundance of H OCO radicals observed in the experiments 
of Oba et al. (2010) suggest that CO2 is readily formed from ex- 
cited HOCO molecules. Therefore, most HOCO radicals con- 
tinue to react before relaxing to their stable state. To reproduce 
these experiments, we assume that 99 % of HOCO radicals are 
sufficiently excited to form CO2, whilst 1% of them are stabi- 
lized; 

IGoumans & Anderssonl (|2010) performed gas-phase O-CO 
potential energy surface calculations, showing that O and CO 
can form a van der Waals complex, allowing O atoms to stay 
bound to CO for a long time. Following these results, we con- 
sider that O atoms that meet CO molecules (via direct accre- 
tion from gas-phase or surface diffusion) form a loosely bound 
O...CO complex. The H atoms that meet these O..CO complexes 
react via a barrierless reaction to form an excited HO. ..CO* com- 
plex. If the time for energy transfer to the surface is long enough, 
the complex can yield OH + CO, or t-HOCO* radical via bar- 
rierless reactions. Otherwise, the complex forms the t-HOCO 
radical through quantum tunneling. We consider that 99% of 
HO. ..CO* complexes continue to react without activation bar- 
riers. 

Appendix B lists all the reactions involved in the formation 
of CO2 with their corresponding activation barriers and trans- 
mission probabilities. 

A.3. Quantum chemical calculations 

The OH + H2 reacti on system has been theoretically studied by 
Nguy en et alJ (1201 ll) . These authors have computed the forward 
and reverse reactions, the imaginary frequency of the transition 
states, and the rate constants of the eight reactions involving H2, 
HD, D2, OH, and OD using the semiclassical transition-state the- 
ory (SCTST). 

Th e H2 O2 + H reaction has b een studied by iKoussa et alJ 
(2006) and Ellings on et al.l (120071) . However, data for the re- 
actions involving deuterated isotopologues was not available. 
Therefore, to obtain the data required for the model, quantum 
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che mistry calculation s have been conducted with the Gaussian 
09(Frisch etal.2009) program. For these calculations, the PBEO 
(Perdew et al. 1996b a; Adamo & Barone fl999h functional and 
the aue-cc-pVTZ (iDunning Jrlll989t iKendall et al] 1 19921) basis 
set were used because this combination produced results that are 
in good agreement with the experimentally det ermined values 
for the process in gas-phase (Klemm e t al.ll 19751 who measured 
an activation barrier of 4.6 kcal/mol = 2300 K) . 

The CO + H reaction has been theoretically studied by sev- 
eral authors (IWoonl 120021 lAndersson etail 1201 It " 



Peter s et aD 



120121) . We decided to use the work of iPeters etal.l (120121) be- 
cause they used the most accurate methodology and obtained a 
value for the activation energy for the formation of HCO that 
best agrees the gas-phase experiment dWang et"al]|1973l) . The 
transmission probabilities of all the reactions producing deuter- 
ated formaldehyde and methanol ar e computed from this reac- 
tion and rela t ive ra tes measured by Naga oka et al.l (120071) and 
iHidaka et al] J2009b or deduced by TCK12b. 

The formation of carbon dioxide includes two rea ctions hav- 
ing an activation barrie r. Quan tum cal culations of | Talbi et al 



(2006&; iGoumans et alj d2008l) . and iGoumans & Anderssoi 
(l2010h showed that reaction (fTTT i has an activation barrier 
of 2500 - 3000 K, leading to a transmission probabili ty o f 
5 x 10~ 23 (IGoumans & Anderssonl2010tlGarrod & Paulvll201 lb . 
Assuming an activation energy of 2500 K and considering a 
square barrier, we reproduce this transmission probability with a 
barrier width of 0.8 A. The transmission probabilities of the re- 
action pathways involved in reaction (fT2l . and including HOCO 
radicals, the van der Waals complex HO... CO, and their deuter- 
ated isotopologue s, were deduced from the potential energy sur- 
face computed bv lYu et aD d200ll) . These authors computed the 
stationary points of the potential energy surface of this reac- 
tion using an extrapolated full coupled cluster/complete basis set 
(FCC/CBS) method. 
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Table B.l: List of grain surface chemical reactions considered in this work along with their activation barriers and transmission 
probabilities. 



Reaction 














E a (K) 


P r 


H 


+ 


H 


— > 


(o)H 2 









0.75 


H 


+ 


H 


-> 


(P)H 2 









0.25 


H 


+ 


D 


-> 


HD 









1 


D 


+ 


D 


— > 


(o)D 2 









0.66 


D 


+ 


D 


— > 


(P)D 2 









0.33 


CO 


+ 


H 


-» 


HCO 






1763 


1.92(-07) 


CO 


+ 


D 


— > 


DCO 






/ 


1.92(-08) 


HCO 


+ 


H 


— > 


H 2 CO 









1 


HCO 


+ 


D 


— > 


HDCO 









1 


DCO 


+ 


H 


-> 


HDCO 









1 


DCO 


+ 


D 


-> 


D 2 CO 









1 


H 2 CO 


+ 


H 


— > 


CH3O 






/ 


9.60(-08) 


H 2 CO 


+ 


D 


— > 


CH 2 DO 






/ 


9.60(-09) 


H 2 CO 


+ 


D 


— > 


HCO 


+ 


HD 


/ 


9.31(-08) 


H 2 CO 


+ 


D 


-> 


HDCO 


+ 


H 


/ 


9.31 (-08) 


HDCO 


+ 


H 


-» 


CH 2 DO 






/ 


1.11 (-07) 


HDCO 


+ 


D 


— > 


CHD 2 






/ 


1.11 (-07) 


HDCO 


+ 


H 


— > 


HCO 


+ 


HD 


/ 


1.54(-07) 


HDCO 


+ 


D 


— > 


DCO 


+ 


HD 


/ 


9.31 (-08) 


HDCO 


+ 


D 


-> 


D 2 CO 


+ 


H 


/ 


9.31(-08) 


D 2 CO 


+ 


H 


-> 


CHD 2 






/ 


1.27(-07) 


D 2 CO 


+ 


D 


— > 


CD3O 






/ 


1.27(-08) 


D 2 CO 


+ 


H 


-> 


DCO 


+ 


HD 


/ 


7.30(-08) 


CH 3 


+ 


H 


-> 


CH3OH 









1 


CH3O 


+ 


D 


— > 


CH3OD 









1 


CH 2 OH 


+ 


H 


— > 


CH3OH 









1 


CH 2 OH 


+ 


D 


— > 


CH 2 DOH 









1 


CH 2 OD 


+ 


H 


-> 


CH3OD 









1 


CH 2 OD 


+ 


D 


— > 


CH 2 DOD 









1 


CH 2 DO 


+ 


H 


— > 


CH 2 DOH 









1 


CH 2 DO 


+ 


D 


-> 


CH 2 DOD 









1 


CHDOH 


+ 


H 


— > 


CH 2 DOH 









1 


CHDOH 


+ 


D 


— > 


CHD 2 OH 









1 


CHDOD 


+ 


H 


-> 


CH 2 DOD 









1 


CHDOD 


+ 


D 


— > 


CHD 2 OD 









1 


CHD 2 


+ 


H 


— > 


CHD 2 OH 









1 


CHD 2 


+ 


D 


— > 


CHD 2 OD 









1 


CD 2 OH 


+ 


H 


-> 


CHD 2 OH 






/ 


1 


CD 2 OH 


+ 


D 


-> 


CD3OH 









1 


CD 2 OD 


+ 


H 


— > 


CHD 2 OD 









1 


CD 2 OD 


+ 


D 


— > 


CD3OD 









1 


CD3O 


+ 


H 


— > 


CD3OH 









1 


CD3O 


+ 


D 


-> 


CD3OD 









1 


CH3OH 


+ 


D 


-> 


CH 2 OH 


+ 


HD 


/ 


2.88(-07) 


CH 2 DOH 


+ 


D 


— > 


CHDOH 


+ 


HD 


/ 


1.92(-07) 


CHD 2 OH 


+ 


D 


-» 


CD 2 OH 


+ 


HD 


/ 


1.50(-07) 


CH3OD 


+ 


D 


-> 


CH 2 OD 


+ 


HD 


/ 


2.88(-07) 


CH 2 DOD 


+ 


D 


-» 


CHDOD 


+ 


HD 


/ 


1.92(-07) 


CHD 2 OD 


+ 


D 


— > 


CD 2 OD 


+ 


HD 


/ 


1.50(-07) 


O 


+ 





— > 


o 2 









1 


o 2 


+ 





-» 


O3 









1 


O 


+ 


H 


-> 


OH 









1 





+ 


D 


— > 


OD 











OH 


+ 


H 


— > 


H 2 











OH 


+ 


D 




HDO 











OD 


+ 


H 




HDO 











OD 


+ 


D 




D 2 











OH 


+ 


OH 


— > 


H 2 2 






(R= 0.8) 




OD 


+ 


OH 




HD0 2 






(R= 0.8) 




OD 


+ 


OD 




D 2 2 






(R= 0.8) 





22 



V. Taquet et al.: Water ice deuteration: a tracer of the chemical history of protostars 



OH 


+ 


OH 




H 2 


+ 





(R=0.2) 






OD 


+ 


OH 




HDO 


+ 





(R=0.2) 






OD 


+ 


OD 


— > 


D 2 


+ 





(R=0.2) 






OH 


+ 


H 2 




H 2 


+ 


H 


2935 


4.07( 


-07) 


OD 


+ 


H 2 




HDO 


+ 


H 


2855 


3.62( 


-07) 


OD 


+ 


HD 




D 2 


+ 


H 


3051 


1.00< 


-09) 


OD 


+ 


HD 




HDO 


+ 


D 


3026 


8.07< 


-10) 


OH 


+ 


HD 


— > 


H 2 


+ 


D 


2789 


8.74( 


-07) 


OH 


+ 


HD 


_> 


HDO 


+ 


H 


2900 


2.8K 


-09) 


OH 


+ 


D 2 


_> 


HDO 


+ 


D 


2703 


7.99( 


-07) 


OD 


+ 


D 2 


_> 


D 2 


+ 


D 


2870 


2.26( 


-09) 


2 


+ 


H 


— > 


H0 2 













O2 


+ 


D 


— > 


D0 2 













H0 2 


+ 


H 




H 2 2 













H0 2 


+ 


D 




HD0 2 













D0 2 


+ 


H 


— > 


HD0 2 













D0 2 


+ 


D 


— > 


D 2 2 













H 2 2 


+ 


H 




H 2 


+ 


OH 


2508 


1.37( 


-07) 


H 2 2 


+ 


D 


— > 


HDO 


+ 


OH 


2355 


5.54( 


-09) 


HDO2 


+ 


H 


— > 


HDO 


+ 


OH 


2523 


1.23( 


-07) 


HDO2 


+ 


H 


— » 


H z O 


+ 


OD 


2524 


1.22( 


-07) 


HDO2 


+ 


D 




D 2 


+ 


OH 


2369 


5.28( 


-09) 


HDO2 


+ 


D 




HDO 


+ 


OD 


2367 


5.29( 


-09) 


D 2 2 


+ 


H 


— > 


D 2 


+ 


OH 


2540 


1.08( 


-07) 


D 2 2 


+ 


D 


— > 


D z O 


+ 


OD 


2384 


4.28( 


-09) 


O3 


+ 


H 


— > 


2 


+ 


OH 









O3 


+ 


D 




2 


+ 


OD 









N 


+ 


H 




NH 













N 


+ 


D 


— > 


ND 













NH 


+ 


H 




NH 2 













NH 


+ 


D 




NHD 













ND 


+ 


H 


— > 


NHD 













ND 


+ 


D 


— > 


ND 2 













NH 2 


+ 


H 


— > 


NH 3 













NH? 


+ 


D 




NH 2 D 













NHD 


+ 


H 




NH 2 D 













NHD 


+ 


D 


— > 


NHD 2 













ND 2 


+ 


H 


— > 


NHD 2 













ND 2 


+ 


D 


— > 


ND 3 













C 


+ 


H 




CH 













C 


+ 


D 




CD 













CH 


+ 


H 


— > 


CH 2 













CH 


+ 


D 




CHD 













CD 


+ 


H 




CHD 













CD 


+ 


D 


— > 


CD 2 













CH 2 


+ 


H 


— > 


CH 3 













CH 2 


+ 


D 


— > 


CH 2 D 













CHD 


+ 


H 




CH 2 D 













CHD 


+ 


D 




CHD 2 













CD 2 


+ 


H 


— > 


CHD 2 













CD 2 


+ 


D 


— > 


CD 3 













CH 3 


+ 


D 


— > 


CH 3 D 













CH 2 D 


+ 


H 




CH 3 D 













CH 2 D 


+ 


D 




CH 2 D 2 













CHD 2 


+ 


H 


— > 


CH 2 D 2 













CHD 2 


+ 


D 


— > 


CHD 3 













CD 3 


+ 


H 




CHD 3 













CD 3 


+ 


D 




CD 4 













HCO 


+ 







C0 2 


+ 


H 









DCO 


+ 





— > 


C0 2 


+ 


D 









CO 


+ 







C0 2 






2500 


4.80( 


-23) 


CO 


+ 


OH 




t-HOCO 






285 


3.50( 


-03) 


CO 


+ 


OD 


— > 


t-DOCO 






121 


1.63( 


-01) 



V. Taquet et al.: Water ice deuteration: a tracer of the chemical history of protostars 



CO 


+ 


OH 




c-HOCO 






2128 


5.81(-17) 


CO 


+ 


OD 




c-DOCO 






1964 


9.57(-16) 


t-HOCO* 


+ 




— > 


CO 


+ 


OH 


13050 


1.63(-10) 


t-HOCO* 


+ 




— > 


c-HOCO 






4114 


1 


t-HOCO 


+ 


H 




co 2 


+ 


H 2 





1 


t-HOCO 


+ 


H 


— > 


H 2 


+ 


CO 





1 


t-HOCO 


+ 


H 


— > 


HCOOH 









1 


t-HOCO 


+ 


D 


— > 


co 2 


+ 


HD 





1 


t-HOCO 


+ 


D 




HDO 


+ 


CO 





1 


t-HOCO 


+ 


D 




HCOOD 









1 


c-HOCO* 


+ 




— > 


t-HOCO 






3272 


1 


c-HOCO* 


+ 




— > 


C0 2 


+ 


H 


12440 


2.63(-01) 


c-HOCO 


+ 






t-HOCO 






3272 


1.19C-17) 


c-HOCO 


+ 






C0 2 


+ 


H 


12440 


2.94(-21) 


c-HOCO 


+ 


H 




C0 2 


+ 


H 2 





1 


c-HOCO 


+ 


H 


— > 


H 2 


+ 


CO 





1 


c-HOCO 


+ 


H 




HCOOH 









1 


c-HOCO 


+ 


D 




C0 2 


+ 


HD 





1 


c-HOCO 


+ 


D 


— > 


HDO 


+ 


CO 





1 


c-HOCO 


+ 


D 


— > 


HCOOD 









1 


t-DOCO* 


+ 




— > 


CO 


+ 


OD 


13220 


9.27(-02) 


t-DOCO* 


+ 






c-DOCO 






4238 


1 


t-DOCO 


+ 


H 




C0 2 


+ 


HD 


(R=0.33) 


1 


t-DOCO 


+ 


H 


— > 


HDO 


+ 


CO 


(R=0.33) 


1 


t-DOCO 


+ 


H 


— > 


DCOOH 






(R=0.33) 


1 


t-DOCO 


+ 


D 


— > 


C0 2 


+ 


D 2 


(R=0.33) 


1 


t-DOCO 


+ 


D 




D 2 


+ 


CO 


(R=0.33) 


1 


t-DOCO 


+ 


D 




DCOOD 






(R=0.33) 


1 


c-DOCO* 


+ 




— > 


t-DOCO 






3343 


1 


c-DOCO* 


+ 




— > 


C0 2 


+ 


D 


13230 


6.18(-02) 


c-DOCO 


+ 






t-DOCO 






3343 


1.51 (-22) 


c-DOCO 


+ 


D 




C0 2 


+ 


D 2 


(R=0.33) 


1 


c-DOCO 


+ 


D 




D 2 


+ 


CO 


(R=0.33) 


1 


c-DOCO 


+ 


D 


— > 


DCOOD 






(R=0.33) 


1 


c-DOCO 


+ 


H 




C0 2 


+ 


HD 


(R=0.33) 


1 


c-DOCO 


+ 


H 


_^ 


HDO 


+ 


CO 


(R=0.33) 


1 


c-DOCO 


+ 


H 


_^ 


DCOOH 






(R=0.33) 


1 


o 


+ 


CO 


_^ 


O...CO 









1 


O...CO 


+ 


H 


— > 


HO.. .CO 









1 


HO.. .CO 








t-HOCO 






775 


9.46(-09) 


HO.. .CO 








c-HOCO 






2618 


8.17(-21) 


HO.. .CO* 






— > 


t-HOCO 






775 


1 


HO.. .CO* 






— > 


c-HOCO 






2618 


1 


DO.. .CO 






— > 


t-DOCO 






727 


2.97(-18) 


DO...CO 








c-DOCO 






2769 


1.60(-20) 


DO...CO* 








t-DOCO 






727 


1 


DO.. .CO* 








c-DOCO 






2769 


1 



